Loop quantization of spherically symmetric midisuperspaces and loop quantum 
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Abstract 

We elaborate on the Ashtekar's formalism for spherically symmetric midisuperspaces and, for loop quan- 
tization, propound a new quantization scheme which yields a graph-preserving Hamiltonian constraint op- 
erator and by which one can impose the fundamental discreteness of loop quantum gravity a la the strategy 
of "improved" dynamics in loop quantum cosmology (LQC). Remarkable consequences are inferred at the 
heuristic level of effective dynamics: first, consistency of the constraint algebra regarding the Hamiltonian 
and diffeomorphism constraints fixes the improved quantization scheme to be of the form reminiscent of the 
improved scheme in LQC which preserves scaling invariance; second, consistency regarding two Hamiltonian 
ly-j constraints further demands the inclusion of higher order holonomy corrections and fixes a ratio factor of 2 

for the improved scheme. It is suggested that the classical singularity is resolved and replaced by a quantum 
bounce which bridges a classical solution to another classical phase. However, the constraints violate briefly 
q during the bouncing period, indicating that one cannot make sense of symmetry reduction by separating 

13^ the degrees of freedom of the full theory into spherical and non-spherical ones in the vicinity of the quan- 

tum bounce, although the heuristic effective dynamics can still give a reliable semiclassical description of 
large-scale physics. Particularly, for the Schwarzschild solution in accordance with the Kruskal coordinates, 
revealing insights lead us to conjecture the complete quantum extension of the Schwarzschild spacetime: 
the black hole is evaporated via the Hawking radiation and meanwhile the quantum spacetime is largely 
extended from the classical one via the quantum bounce, suggesting that the information paradox might be 
resolved. 
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I. INTRODUCTION 



Loop quantum gravity (LQG) [IH3] is considered to be a promising candidate for quantum 
theory of gravity. It provides mathematically rigorous foundation for quantum gravity and has led 
to several significant results about the quantum structure of space and time. While its kinematics 
is well understood, many open issues remain unsettled, mainly in regard to the dynamics and 
the semiclassical limit. In spite of difficulties in the full theory of LQG, the loop approach has 
been successfully applied for symmetry-reduced minisuperspaces in the context of loop quantum 
cosmology (LQC). LQC has yielded valuable insight into the full theory as many issues are better 
understood in LQC thanks to its formal simplicity. In recent years, LQC has undergone lively 
progress and become an active area of research [3H6] . Markedly, it has been shown that, for a variety 
of models of LQC, the cosmological singularity (big bang, big crunch, big rip [7j, etc.) is resolved 
and replaced by the big bounce, therefore affirming the long-held conviction that singularities in 
general relativity signal a breakdown of the classical theory and should be resolved by the quantum 
effects of gravity. 

It is natural to ask whether the black hole singularity is resolved as well. To study loop 
quantum geometry of black holes, the simplest step is to consider the interior of a Schwarzschild 
black hole, in which the metric components are homogeneous with Kantowski-Sachs symmetry. 
By virtue of homogeneity, the loop quantization of the Schwarzschild interior can be formulated 
as a minisuperspace model in a similar fashion to LQC. This has been developed in [BrllOj and 
its effective solution has been investigated at the level of heuristic dynamics in based on the 
original quantization strategy (also referred to as "//o-scheme" ) and in [121CE3] based on the modified 
quantization strategies ("/2-scheme" and "/t'-scheme"). In the hq- and /2-schemes, the black hole 
singularity is resolved by the quantum bounce, which bridges the black hole interior with a white 
hole interior. In the /i'-scheme, the black hole singularity is resolved and the event horizon is also 
diffused by the quantum bounce, through which the black hole is connected to a baby black hole 
with a much smaller mass. 

Both jiQ- and //-schemes suffer from the problem that the resulting dynamics depends on the 
arbitrary choice of the finite sized comoving cell to which the spatial integration is restricted to 
make the Hamiltonian finite. On the other hand, the //'-scheme is independent of the choice of 
the finite sized cell and thus considered to be the correct quantization strategy in accordance 
with the idea of "improved dynamics" first suggested in |14j and later generalized for Bianchi I 
models in |15l I16j . As opposed to the results of the //'-scheme heuristic dynamics, however, the 
spacetime curvature in the immediate vicinity of the event horizon can be fairly flat and does not 
necessarily incur any quantum corrections. The reason that the event horizon wrongly receives 
quantum corrections in the //'-scheme is essentially because the finite sized cell collapses as its side 
surface shrinks when approaching the horizon. The collapse near the horizon is merely an artifact 
resulting from coordinate singularity; in actual fact, a given finite sized cell will simply pass through 
the horizon without being collapsed. This glitch calls into question the minisuperspace treatment 
for the Schwarzschild interior, even though it is legitimate for the Kantowski-Sachs cosmology as 
studied in [17]. 1 Therefore, one is obliged to consider the black hole interior and exterior as a whole 
in the framework of spherically symmetric midisuperspaces. 

Loop theories of spherically symmetric midisuperspaces are of theoretical importance in its 
own right, as they deal with the simplest field-theory framework which is symmetry-reduced. They 



1 The Schwarzschild interior cannot be considered as a self-contained university, as it is extensible beyond its bound- 
ary (event horizon). Put differently, classical general relativity does not admit a vacuum cosmological solution with 
Kantowski-Sachs symmetry. With inclusion of matter content, on the other hand, the Kantowski-Sachs cosmology 
exists and the LQC treatment based on it is sensible (see [17] for the case with a scalar field). 
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provide an arena for testing important issues which are too difficult in the full theory and trivialized 
in minisuperspace (LQC) models. Particularly, the SU(2) internal gauge is reduced to U(l) and 
the 3-dimensional diffeomorphism is reduced to 1-dimensional, making the constraint algebra much 
simpler yet nontrivial. Ashtekar's formalism for spherically symmetric midisuperspaces and its loop 
quantization have been studied and developed with different degrees of rigor [T8H2T] . In this paper, 
we recapitulate results of previous works and elaborate on the geometrical meaning of the Ashtekar 
variables and constraint algebra. For loop quantization, we suggest that some details have to be 
modified and propound a new quantization scheme which yields a graph-preserving Hamiltonian 
constraint operator and by which one can impose the fundamental discreteness of LQG by hand d 
la the strategy of improved dynamics in LQC. 

Rigorous and complete construction for the loop quantum theory of spherically symmetric 
midisuperspaces is faced with complications and still challenging. Nevertheless, ramifications of 
loop quantization and other open issues can still be inferred at the heuristic level of semiclassical 
dynamics. Particularly, an effective solution corresponding to the semiclassical theory has been 
obtained for the complete Schwarzschild spacetime which covers both the interior and exterior 
in [22], yielding a singularity- free global structure akin to that in the /iq- and /2-schemes of the 
minisuperspace (interior) treatment in |llH13j . 

The formulation of [22] follows the strategy of |21j to partially fix the 1-dimensional diffeomor- 
phism gauge in such a way that one is left with a single abelian constraint and a true Hamiltonian. 
The partial gauge fixing avoids the hard problem of having structure functions in the constraint 
algebra and thus makes the consistent loop quantum theory possible, but it also obscures signif- 
icance of the interplay between the diffeomorphism and Hamiltonian constraints. 2 In this paper, 
we do not fix the diffeomorphism gauge but instead propose a new quantization scheme. The new 
quantization scheme does not resolve the complications of constructing a consistent loop quantum 
theory but nevertheless it yields profound insights at the heuristic level of semiclassical dynamics. 
Remarkably, consistency of constraint algebra in heuristic effective dynamics leads to fascinating 
consequences. First, requirement that the Poisson bracket between the Hamiltonian and diffeomor- 
phism constraints has to weakly vanish fixes the improved quantization scheme to be of the same 
form as the //'-scheme in minisuperspace (LQC) models |13 ^ [T5l \TT\ . Second, requirement that the 
Poisson bracket between any two Hamiltonian constraints (with different lapse functions) also has 
to weakly vanish further demands the inclusion of higher order holonomy corrections, which was 
first suggested in [23] [25] in the context of LQC but motivated differently. The second requirement 
also fixes a ratio factor of 2 for the improved quantization scheme, which remains ambiguous (and 
is usually set to be 1) in homogenous models. 

Analysis of the heuristic effective dynamics suggests that, with loop quantum corrections, the 
black hole singularity is resolved and replaced by a quantum bounce. The quantum bounce is 
expected to bridge a classical solution to another classical one, as the constraint algebra is exactly 
satisfied before and after the bounce, provided that the improved quantization scheme is adopted 
and higher order holonomy corrections are included to the order of infinity. The constraint algebra, 
however, breaks down briefly during the transition period of the bounce, implying that the space- 
time is no longer invariant under spatial diffeomorphism and change of spacetime foliation after 
the bounce. Nevertheless, the violation is minuscule at large scale and thus the heuristic effective 
dynamics remains reliable after the bounce insofar as large-scale physics is concerned. During the 
transition of the quantum bounce, the intimate matching between diffeomorphism and Hamiltonian 



2 Analogously, on top of the strategy of [22], if one further fixes the diffeomorphism gauge completely before quan- 
tization, one is led to the quantization of Kuchaf [23], which does not offer insights about the structure of the 
quantum spacetime and the question of singularities. 
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constraints is spoiled as a consequence of imposing fundamental discreteness. 3 This indicates that, 
in the full theory of LQG, a given weave state (coherent state of spin networks) which represents 
a smooth, spherically symmetric space at large scale will eventually manifest granularity of spin 
networks when approaching the quantum bounce. In other words, close to the quantum bounce, 
the weave state becomes very granular to the extent that one can no longer make sense of symmetry 
reduction by separating out spherical degrees from non-spherical fluctuations. After the bounce, 
however, the weave state evolves to become smooth and spherically symmetric again. 

To study the interior and exterior of the Schwarzschild black hole at the same time, we choose 
the Kruskal coordinates, which cover the entire manifold of the maximally extended Schwarzschild 
spacetime (i.e. Kruskal extension of the Schwarzschild spacetime) |26|, I27j. Although the numer- 
ical method remains extremely challenging, inspection of the Hamilton's equation at the level of 
heuristic effective dynamics still offers revealing insight about the characteristic natures of the loop 
quantum geometry of the Schwarzschild spacetime. Firstly, the classical singularity is resolved and 
replaced by a quantum bounce which bridges the black/white hole interior to a different classical 
phase. Secondly, inclusion of the Hawking radiation is mandatory and thus the Hawking evapora- 
tion is expected. These lead us to conjecture that the complete quantum Schwarzschild spacetime 
manifests the similar structure of the 2-dimensional quantum black holes studied in [28-30J: the 
black hole is evaporated via the Hawking radiation, and the quantum spacetime is largely extended 
from the classical one as the black hole interior is extended beyond the putative singularity. 4 The 
information paradox due to the Hawking evaporation is believed to be resolved by the quantum 
extension of the classical spacetime as for the 2-dimensional black holes in |28H30| . 

This paper is organized as follows. In Sec. [ll| we set up Ashtekar's formalism for spherically sym- 
metric spacetimes. In Sec. Ill, the classical solution in terms Ashtekar variables for the maximally 
extended Schwarzschild spacetime is given in accordance with the Kruskal coordinates. In Sec. |IV[ 
loop quantization of spherically symmetric midisuperspaces is discussed. In Sec. [Vj we study the 
heuristic effective dynamics and conjecture the quantum geometry of the Schwarzschild spacetime. 
We end with a summary and discussion in Sec. |VI[ Supplementary materials are included in the 
appendices. 



II. ASHTEKAR'S FORMALISM FOR SPHERICALLY SYMMETRIC SPACETIMES 

In this section, we formulate the canonical theory of gravity in terms of Ashtekar variables 
for the spacetimes which admit spherical symmetry. We mainly follow the previous works in 
|X8|, [20l [21] but give more details and expound on the physical meaning of the canonical variables 
and constraints. 



A. Ashtekar variables and constraints 



The ADM formalism supposes that the 4-dimensional spacetime M is foliated into a family of 
spacelike surfaces T, t , labeled by the parameter time coordinate t. Further, we consider the cases 



By contrast, in the context of LQC, imposing discreteness by hand yields no inconsistency, because diffeomorphism 
is removed by homogeneity and one is left with scaling invariance (with respect to the choice of the finite sized 
cell). However, as a consequence of imposing discreteness, scaling invariance is softly broken in the quantum theory 
(the fundamental discreteness enters as a constant step size in the difference equation of Hamiltonian constraint, 
while triad variables scale up with the finite sized cell). 

The global structure of the quantum Schwarzschild spacetime conjectured in this paper is very different from that 
obtained in 22 . The discrepancy will be discussed in the end of Sec. VC 
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that the topology of each of the "leaves" is of the form £ = I x S 2 , where I is a 1-dimensional 
manifold coordinatized by x (called "radial coordinate") and S 2 is a 2-sphere coordinatized by 9 
and (p. 5 

If E admits spherical symmetry (i.e., there are three Killing vectors of S 2 that form the so(3) = 
su(2) algebra), the Ashtekar connection in the full theory is symmetry-reduced to the form of 

A = A i r i dx a = A x (x)t3<1x + (Ai(x)t 1 + A 2 {x)T 2 )dO 

+ [(A 1 (x)t 2 - A 2 {x)n) sin 9 + t 3 cos 6]dc/), (2.1) 

where A x , A± and A 2 are functions of x and Tj and = — zcij/2 with <7j being the Pauli matrices are 
orthonormal generators of SU{2) (note n A Tj = fa, 73] = €jj fc Tfc). Correspondingly, the densitized 
triad takes the form 

E = E«T i d a = E x (x)T 3 sm9— + (E\x)T 1 +E 2 (x)T 2 )sm9— 

+ {E\x)t 2 -E 2 {x)t 1 )^ (2.2) 

where again E x , E 1 and E 2 are functions of x. Note that the full (3+1 dimensional) theory has 
been reduced by the spherical symmetry to a 1+1 dimensional theory. In the full theory, the 
connection is a one- form and the densitized triad is a vector density of weight 1. But in the one 
dimension, under the coordinate transformation x —> x(x), A x transforms as a scalar density of 
weight 1, i.e. A x (x) = (dx/dx)A x (x), while A\ and A 2 are scalars; and E x transforms as a scalar, 
while E 1 and E 2 are scalar densities of weight l. 6 

In terms of the symmetry-reduced variables, the symplectic structure given by the symplectic 
two- form reads as: 

= I dx I d9d(f)s\n9 (5E X A 5A X + 2 (5E 1 A 5A X + 8E 2 A SA 2 )) 

8ttGj Jj J 

= — ^- f dx (5E X A SA X + 2 (5.E 1 A <5^i + A <L4 2 )) , (2.4) 

2Cr7 J J 



which implies 



{A x (x),E x (x')} = 2G-f5(x-x'), (2.5a) 
{A 1 {x),E 1 (x')} = G~f5{x-x'), (2.5b) 
{A 2 {x),E 2 (x')} = G~f5(x-x'), (2.5c) 



This covers a wide range of important cases. For the exterior of the Schwarzschild solution, we have x G [0, oo) and 
x — corresponds to the horizon [21) - For the interior of the Schwarzschild black hole [13] or for the Kantowski- 
Sachs cosmology [171. we have x G (— oo, oo) and one can further reduce the phase space to that of finitely many 
degrees of freedom by imposing homogeneity. In this paper, we take into account both the interior and exterior of 
the Schwarzschild solution as a whole, and choose x G (—00,00) to be the horizontal coordinates in the Kruskal 
diagram. 

It is helpful to keep track of the density character of the canonical variables. For example, it can be used to check 



that U(x) defined in ( |2.10[ ), D(x) in ( |2.12[ ) and C{x 

1, respectively. (Correspondingly, X(x), N x (x) and N(x) are a scalar, a scalar density of weight 



in (2.151 are scaler densities of weight 1, weight 2, and weight 

and a scalar, 



respectively). Furthermore, The minus sign in the integrand in ( |2.12[ ) is due to the fact that A x is a scalar density 
and E x is a scalar, while the pairs (Ai,^ 1 ) and (A 2 ,E 2 ) are the other way around. Also see Footnote [9] and 
Appendix |A] 
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and all the other brackets vanish. G is Newton's constant and 7 is the Barbero-Immirzi parameter. 



The connection in (2.1) gives the field strength 7 

F = ^F ab dx a Adx b = dA + AAA 

= {A\ + A\ - l) r 3 sin 9 d6 A d<f> - [( A[t 2 - A 2 n ) -A x (Am + A 2 t 2 )} sin 9d(j)Adx 

+ [( A[n + A' 2 t 2 )+A x (A lT2 - A 2 Ti)] dx A d9. (2.6) 

In the full theory [I], the dynamics is described by three constraints. First, the Gauss constraint 
function is given by 

C G [K\x)] := — J tPxVGu (2.7) 

which generates the internal SU(2) gauge transformation with 

d = V a E? := d a E? + e l3 k A{E a k (2.8) 



for any smooth field A(x) on E. With the symmetry-reduced variables in (2.1) and (2.2), it is easy 
to show that the two components Gi=\ and G{ =2 vanish identically and the only non-vanishing 
component is 

G i=3 = sin 9 [E xl + 2AXE 2 - 2A 2 E l ) . (2.9) 
Thus, the SU(2) gauge is reduced to U(l) and generated by 

Cg[A(x)] = [ dxX (E x/ + 2A 1 E 2 - 2A 2 E l ) =: / dxX(x)U(x), (2.10) 

2<^7 J I J 1 

where A(x) := xr J d9d(j) sin 9 A 7,=3 (x, 9, <f>) and note that the factor (47r) _1 is introduced to coun- 
terbalance the angular integral. 8 



7 Throughout this paper, we will use ' and ' to denote derivatives with respect to the radial coordinate x and the 
parameter time t, respectively. However, in some occasions, ' is merely used to denote a related but different 
variable, such as S(x — x') as x is in contrast to x or E — > E' as E' is transformed from E. The meaning of ' 
should have no confusion when put in the context. 



Under the 17(1) gauge transformation, the term 2A\E 2 — 2A2E 1 in (2.101 corresponds to a rotation on the (ri,T2 



plane, while E x gives rise to the inhomogeneous term for the transformation of A. That is 
5E\x) = ^(x),^ ^dx'\(x')(2A 1 (x')E 2 (x')-2A 2 (x')E 1 (x'))^ 

and similarly 

8E 2 (x) = \{x)E 1 {x), 5Ai{x) = -X(x)A 2 (x), SA 2 (x) = X(x)Ai(x), 

while 

5Ax(x) = 
and obviously 

5E x (x)=0. 

To summarize, the SU(2) gauge transformation e T ' A( - x ^ is reduced to the U (1) transformation e T3A< - :E ' , which rotates 
(tlj T2) by the angle X(x), and the fields transform accordingly: 

E^fE 1 = e T3X(x) Ee- T 3A(x) 

and 

A -> A' = e T3Hx) Ae- T3K ' x) + e T3Kx) de~ T3X{x) = e ^^) Ae -^K^) _ dx \( x ) Tzdx . 



{^^'2^7 f dx '^( x '^x' EX ( x '^ = J dx'X(x')d x i&(x — x) = -d x X(x) 
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Second, the diffeomorphism constraint function is given by 

C Di ft[N a (x)} := g-L- J d 3 x (N a E?Fi b - N a A\G t ) , (2.11) 

which generates diffeomorphisms in £ for any smooth vector field N a (x) . In terms of the symmetry- 
reduced variables, again, two components of the integrand coupled with N e and vanish iden- 
tically and the only non-vanishing contribution is given by 

C ms [N x (x)} = -L f dxN x (2A\E l + 2A' 2 E 2 - A X E X ') =: [ dxN x {x)D(x), (2.12) 
2<j7 J j J j 

where N x (x) := / d8d<ft sin 9 N a=x (x, 9, (/)). This generates the remnant diffeomorphism along 

Q 

x. 

Finally, the Hamiltonian constraint (also called scalar constraint) is given by 

° [Nm := itc L d3xN ^ {^ k K b E]E b k - 2(1 + ^KlK^EfEf] , (2.13) 

I ~ 1 1/2 

where /[ a / 6 ] := (f a f b - /&/ a )/2 and e := |detEf| (note that det Ef = detq ab =: q by the 



relation qq ab = EfE b -5^). The triad in (2.2) gives detEf = E X [(E 1 ) 2 + (E 2 ) 2 ] sin 2 9. The extrinsic 
curvature takes the form 

K = KiTidx a = K x (x)T 3 dx + {K x (x)ti + K 2 (x)r 2 ) d9 + (Ki(x)r 2 - K 2 (x)n) sin0#, (2.14) 



9 Under an infinitesimal diffeomorphism transformation x — > x = x + e(x), A\ transforms as a scalar: Ai(x) 
A^x) = Ai(x). Thus, 

6Ai(x) = Ai(x) - A(x) = Ai(x - e(x)) - A x {x) = Ai(x) - e(x)d$Ai(x) - Ai(ar) + C(e 2 ) 
= -c(ar)fl! r i4i(!c) + C(e 2 )- 

And E 1 transforms as a scalar density: E 1 (x) -> = (dx/dij.E 1 (as) = (1 - ft,e(a:)) E 1 (x) + 0(e 2 ). Thus, 

5E X (x) = E 1 (x)-E 1 (x) = E 1 (x)-e(x)d x E 1 (x)-E 1 (x) + 0(e 2 ) 
= -9 x e(a;)E 1 (a:) - e(x)d x E 1 (x) + 0(e 2 ). 



This is in agreement with the 2A'^E term in (2.12 l provided that x — x — N x (x)dt, as 

SAiQc) = $Ai(x), ^ £dx'2N x (x)A , 1 (x')E 1 (x')j = J dx N x (x')Ai{x')8{x - x) = N x (x)d x Ai(x), 

8E 1 (x) = i^ l {x),^J^dx'2N x {x')A' 1 {x')E 1 {x')\=-J dx' N x {x')d x ,8{x - x')E\x) 

= j dx'8(x - x')d x , (N x {x')E 1 {x)) = d x N x {x)E 1 {x) + N x {x)d x E 1 {x) 

Similarly, 2A' 2 E 2 term corresponds to the diffeomorphism transformation for the pair (A2, E 2 ). On the other hand, 
because (A X ,E X ) has opposite density character, it corresponds to —A x E xl with a minus sign (recall Footnote Rjb. 
Generically, for any scalar density f(x) of density weight w, i.e. f(x) — > f(x) — (dx / dx) w f (x) , we have 

8f(x) = -wd x e{x)f(x) - e{x)d x f{x) + 0{e 2 ), 

which is in agreement with 

8f(x) = {f(x),C DiB [N x ]} = wd x N x (x)f{x) + N x {x)d x f{x) =: £*»/(*), 
if /(a;) is a composite function of the canonical variables. 
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where K x , K\ and K 2 can be written in terms of A's and E f s. By (2.2) and (2.6), the reduced 
variables then yield 



C[N(x)} 



1 



dx- 



N(x) 



2Gjj y/i&wW^TWF] 

x^2E x {E 1 A' 2 - E 2 A[) + 2A X E X (A^ 1 + A 2 E 2 ) + {A\ + A\ - l)^ 1 ) 2 + (E 2 ) 2 } 
-(1 + 7 2 ) (2K X E X (K 1 E 1 + K 2 E 2 ) + [{K^ 2 + {K 2 ) 2 ][(E 1 ) 2 + (E 2 ) 2 }) } 

dxN(x)C(x), (2.15) 



where JV(x) = 4^/ <Z0# sin 9 N(x, 9, 



B. Polar type variables 

In what follows, to be better adapted to the U(l) gauge transformation, we follow the ideas in 
[THEmE] to introduce the "polar" type variables. 10 First, we define A p and E p via 



A\ = A p cos (3, A 2 = Apsin/3, 

E 1 = E p cos{a + p), E 2 = E p sin(a + (3). 



(2.16a) 
(2.16b) 







Note that ^4 p , -E' p and the inner product (Ai,A 2 ) • (E 1 ^ 2 ) = A p E p cos a are invariant under the 
U(l) rotation; i.e. the angle a(x) between (Ai,A 2 ) and (E 1 ^ 2 ) is gauge independent, while the 
angle (3(x) is pure gauge. With the new variables, (2.3) becomes 

— - / dx\5E x A SA X + 2( cos a 5E P A 5A p + sin a A p 5E p A 5f3 

+ sin a E p 5A p A 5(3 -\- sin a E p 5A p A <5a + cos a A p E p 5a A 5/3) } 

— ^- / dzj&E 1 A (54,. + 2 [$(£'' cos a) A 5 A p + 5 {A p E p sin a) A5p]\ (2.17a) 

— ^- [ dx\5E x A5A X + 2[5E P A 5(A p cos a) + 5{A p E p sin a) A5(a + (3)]\. (2.17b) 



Equation (2.17b) tells that one can transform the set of canonical variables (A x , A\,A% E x , E 1 , E 2 ) 
to the new one (A x , A p , 77; E x , E p , F" 1 ) by defining 



n 



2A p cos a, 



a + (3, 

2A p E p sina = 2A 1 E 2 - 2A 2 E X . 



Equivalently, we have 



{A x (x),E x (x')} = 2G 7 S(x-x'), 
{A p (x),E p (x')} = 2Gj5(x-x'), 
{rj(x),P r '(x')} = 2Gj8(x-x'), 



(2.18a) 
(2.18b) 
(2.18c) 



(2.19a) 
(2.19b) 
(2.19c) 



10 In this paper, we adopt slightly different notations for polar type variables, which we believe are better than those 
used in [TM22] . See Appendix 111 
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and all other brackets for the new canonical variables vanish. 11 



In terms of the polar type variables, the Gauss constraint (2.10) reads as 
1 



U 



2G 7 



E xl + 2A X E 2 - 2A 2 E X ) = ^- {E xl + 2A p E" sin a) = (E xl + P r >) , (2.20) 



and the diffeomorphism constraint (|2.12|) becomes 
D 



1 



2A\E X + 2A 2 E' Z - A X E 



I 77.2 



2G 7 
1 

1 

2Chf 



1 



2G 7 



2G 7 



(2A' p E p cos a + 2A p E p sin a/3' - A X E X ') 



[((2A p cosa)' + 2A p sin a c/) E p + 2A p E p sin a/3' — A x E a 
(A' p E p + j//™ - . 



(2.21) 



It is more involved to write the Hamiltonian constraint in terms of canonical variables as we 
first have to solve (K x , Ki, K 2 ) in terms of A's and E's. To make the calculation tractable, we 



perform change of variables in accordance with (2.16b) by rotating the orthonormal basis (tl, t 2 , t 3 ) 



of su(2) to the new orthonormal basis {t\i ,t 2 i , T3) (orientation unchanged): 

Ti'(x) := ti cos(q + f3) + t 2 sin(a + /?), 
t 2 i{x) := — T\ sin(a + (3) + T2 cos(a + /?)■ 



(2.22a) 
(2.22b) 



It should be noted that the rotation angle is different from point to point in x. Both t\> and t 2 i 
are functions of x and we have 



t[, = t 2 , (a + f3)' , t' 2 , = -T V (a + P)' . 



(2.23) 



Equations (2.1) and (2.2) then read as 



A = A l a Tidx a = A x (x)T3dx + A p (x) (cos a Ty (x) — sin a t 2 i(x)) d6 

+A p (x) (sin a Ty(x) + cos aT 2 /(x)) sin 9 dcj) + T3 cos < 



(2.24) 



E = EtT { d a = E x (x)T 3 sm0^- + E p (x)T V (x)sm9-^+E p (x)T 2 ,(x)^, (2.25) 

as Ef is diagonalized in the new basis. We have q = detEf = E X (E P ) 2 sin 2 6 and the relation 
e a = ykl (Ef) -1 yields the (undensitized) cotriad: 



e = e l a Tidx a 



E p 



: T 3 dx + y/\E x \ T\i (x)d6 + \f\E*\ sin r 2 / {x)d(f). 



(2.26) 



Correspondingly, the 3-metric components given by q a b = e„e£#y read as 



[EPf 
\E T \ 



q ee = \E x l q H = \E X \ sin 2 1 



(2.27) 



11 Alternatively, by (2. 17a I, one can transform the original canonical variables to a different canonical set 
(A x ,A p ,f3;E x ,P p ,P p ) with P p := 2£; p cos« and := 2A p E p sina (see [THUD])- This canonical transformation 
is less relevant for our purpose. 
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The cotriad in (2.26) gives 



de 



E x 



2y/W\ 

+ 



Tydx A d6 + \/\E*\ (a + P)'T 2 'dx A dd 
T2> sin 6 dx A dip — \J \ E x \ (a + B)'ry sin 9 dx A- 



+ \f\E x \ t 2 > cos 6d9 A 



(2.28) 



cotriad via de + T A e = is then given by the solution: 



where (2.23) has been used to compute dry and dry- The spin connection V compatible with the 

(2.29) 



E x/ E x/ 

T = T l a Tida = -{a + B)'T 3 dx + T 2 *d9 - -^- p Ty sin 9 dip + r 3 cos I 



Equations (2.24) and (2.29) then give the extrinsic curvature = A — T: 



jK = -fK i a T i dx a 
[A x + (a + p)'] r 3 dx 



(2.30) 



+ 



A p cos a t\i — ( sin a + 



£2/ 

2E^ 



72' 



A p cos ary + ( A p sin a + 



2E~p Tl ' 



sm( 



which takes the form of (2.14) with 





= A x + (a + ft)', 




E xJ 


n/Ki 


= A cos 3 H — - 

p H 2EP 




E xl 


lK 2 


= A a sin 8 — - 

p 2EP 



sin(a + B), 
cos(a + B). 



If we immediately impose the Gauss constraint (2.20), we have 



(2.31a) 
(2.31b) 

(2.31c) 

~fK = ^K i a T i dx a = [A x + (a + B)'] r 3 dx + A p cos a T V dO + A p cos a t 2 , sin 6dcp, (2.32) 

and correspondingly yKi = A p cosa cos(a + B) and "yK 2 = A p cosa sm(a + B). 12 

Remark. In the above calculation, the reader might be puzzled why the basis (ti',t 2 i,t 3 ) is 
treated as rc-dependent while the original basis (n, t 2 , t 3 ) is considered fixed with respect to exterior 
derivative. Is not the difference between these two merely a gauge transformation? Well, yes and 
no! What we just did is not a gauge transformation. In fact we have stuck with the (t\,t 2 ,t 3 ) 
basis all the way through; introducing (ri', t 2 ', t 3 ) is only a procedure for change of variables which 
largely simplifies the calculation (mainly, by diagonalizing Ef). Alternatively, if one wish, one can 
indeed regard the change from (ti,t 2 ,t 3 ) to (ry ,t 2 > ,t 3 ) as a gauge transformation rather than a 
change of variables. This way, (ry , t% , t 3 ) is a fixed basis just like (ti,t 2 ,t 3 ) and we should not 



take dry and dr 2 i as we did in (2.28); consequently, we will not have (a + B)' in (2.28) and (2.29). 



However, under the gauge transformation (2.22), all the fields have to transform accordingly: 
i.e., E -> E' = e - T ^ a+ ^Ee T ^ a+ ^ and A -»• A' = e - T ^ a+ ^Ae T3( - a+ ^ + e - Ts( - a+ ^de T3( - a+ ^ (see 
Footnote pi) . Note that the inhomogeneous term for the above gauge transformation of A yields 



(a + 8)'T 3 dx, which is missing in (2.24). Therefore, both methods (change of variables and gauge 



Equations (2.30 1 and (2.321 give rise to the equivalent Hamiltonian constraint modulo the Gauss constraint. In 



this paper, (2.301 is used to derive (2.42c I 
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transformation) yield the same result in the end; particularly, the term (a + (3)'T^dx for = A — T 
is from T in the former approach and from A in the latter. (It should be noted that both the 
Ashtekar connection A and the spin connection T do not transform covariantly under the gauge 
transformation, but their difference jK does!) 

From (2.16) and (2.31), we can compute the part inside the curly parenthesis of (2.15): 

2Gy/qC := 2E X (E 1 A' 2 - E 2 A[) + 2A X E X (AiE 1 + A 2 E 2 ) + (A 2 + A 2 - 1)[(^ 1 ) 2 + (E 2 ) 2 ] 
-(1 + 7 2 ) {2K X E X (K 1 E 1 + K 2 E 2 ) + [{K x f + (K 2 ) 2 ][{E 1 ) 2 + (E 2 ) 2 }) 

= --^ {4 ( 7 2 + ,4 2 ) (E p ) 2 + (1 + 1 2 )(E X ') 2 + 8 7 2 E x E p A' p sin a + 8E X E? A X A P cos a 

+8E x E p A p cos a ((1 + 7 2 )a' + tf) + 4(1 + j 2 )E p A p E xf sin a}. (2.33) 

Although yJqC as a whole is U(l) invariant, the above expression still involves [/(l)-dependent 
variables A x and (3. To make the U(l) invariance explicit, it is convenient to add a suitable multiple 
of the Gauss constraint: 

2G>yE x d x U = 2E x E p A' p sin a + 2E x E p A p cos a a' + 2E x A p E pl sin a + E X E X ", (2.34) 
so that the new y/qC reads as 

2G^C -> 2G^C + 2GjE x d x U 

i I {4( 7 2 + ^ 2 ) (^) 2 + (l + 7 2 )(^') 2 



4 7 2 



+8£ :c J EM p (A c + a' + /?') cos a + 4(1 + 7 2 )£M p £ x/ sin a} 



(2.35) 



in which A x and f3 appear only through the [/(l)-invariant quantity yK x = A x + (a+(3)' . Replacing 

(2.36) 



A p cos a and A p sina with (2.18) and A 2 with 



A 2 

A 2 „ = -r + A 2 sin 2 a 
" 4 " 



pv 



A 2 

4 V 2EP 



we then have 



IGyfqC 



E x E p 'P ri 



-{E p f 



4 7 2 

l + 7 2 
2 7 2 



r 

1 4- -v 2 

Itfc'pn _ —J. (E xl ) 2 



A 2 p {E p ) 2 E x E p A p 



4 7 2 



7 



4 7 2 

, (A, + ol + /?') 



+ E X E X ", 

E x E xl E pl 



EP 



+ K L + e x E x ", 



(2.37a) 



(2.37b) 



wh ere in the second step we have applied the Gauss constraint G = by the form of (2.20). Finally, 
by |2 15| and noting that y/q = ^J\E X \ [(E 1 ) 2 + (E 2 ) 2 } = E p , we have 



C 



EP 



A 2 p EP 



2G 1 y/lE^ 4 7 2 yp^ 



sgn (E x )^^(A x + a' + P') 



, „ \/W\ E X 'EP' (E x ') 2 , VX JW\E X " 

-sgn(£*) VI ' + - S=^- + sgn{E x ) V 



(2.38) 



(EP) 2 ' 4y/\E?\EP ' ' EP j 

which is now written all in terms of U(l) gauge independent variables: E p , A p , E x , and A x +a' . 
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C. Further change of variables 



Previously, we have transformed the set of canonical variables (A x , A\, A2; E x , E , E 2 ) to the 
new set (A x , A p , 77; E X ,E P , P^). To better grasp the physical meaning of Ashtekar variables, we can 
perform change of variables once more to separate ?7(l)-dependent variables from U (l)-independent 
ones. The symplectic form (2.17b) can be recast as 

= / dx\sE x A 5{A X + a + tf) + 5E P A 5{2A p cos a) 

+6(2A p E p sin a + E xl ) A 5(a + 0) } (2.39) 

by integration by parts. This allows one to canonically transform (A x , A„, 77; E x , E p , P^) to the new 
set of canonical variables (A x , A p ,rj; E x , E p , P^) by defining 



A x 



A x + V ' = A x + (a + p)', 

pr, + E xl_ 



In terms of the new canonical variables (A x , A p , 77; E x , E p , P„), we have 

{A x (x),E x (x')} = 2G 1 8(x-x'), 
{A p (x),E p (x')} = 2Gj6(x-x'), 
{r]{x),P r >{x')} = 2Gj5(x-x'), 



(2.40a) 
(2.40b) 



(2.41a) 
(2.41b) 
(2.41c) 



and all other brackets vanish. Meanwhile, the three constraints given by (2.20), (2.21) and (2.38) 
are recast as 



U 
D 

C 



1 



2G 7 
1 

2Chi 
1 f 



A' p E p + 7]'P r ' - A X E X ') , 



E p 



A 2 p E p 



2G 1 y^Wl 4 7 2 ^^ 



sgn(^) 



y/\&\A p A x 



r 



t ^^/W\E x 'E p ' 
-sgn(E x ) v 



(E 



xi\2 



(E p ) 2 ' 4y/\EP\EP 
The Poisson brackets between the constraints are given by 



+ sgn(£ ,a 



ep~ 



{C G [X],C G [X']} 

{C G [X],C{N]} 
{C G [X],C» m [N x ]} 
{C DiS [M%C DiS [N x ]} 

{C[N],CDis[N x ]} 
{C[N],C[M}} 



0, 

0, 

-C G [N X X'} = -C G [£ m X], 
C Dm [M x N x ' - N X M X '} = C DiS [[M x , N 1 
Cr>ig[£M*N x ] = — CDiff[£NxM x ] 
-C[N X N'} = -C[£ N *N], 



CDiff 



E 2 

(NM' - MN')-r^ 



C G 



E 2 

t/\ ^X 



(NM' - MiVO-rrV 

\q\ 



(2.42a) 
(2.42b) 



(2.42c) 



(2.43a) 
(2.43b) 
(2.43c) 

(2.43d) 
(2.43e) 

(2.43f) 
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The right hand sides of (2.43) are weakly zero (abbreviated as ~ 0) in the sense that they vanish 
in the constrained phase space. 13 The three constraints are of first class in Dirac's terminology, 
but the constraint algebra is open in the sense that we have structure function instead of structure 



constants, as the smearing fields in the right hand side of (2.43f) depend on dynamical fields. (See 
[1] for the constraint algebra in the full theory and related comments.) 

The Hamiltonian in the full theory is given by % = C[N(x)\ + Ci)is[N a (x)] + Cg[uj 1 • t(x)] (see 
P]). With the imposition of spherical symmetry, T~L reads as 



U = C[N]+C DiS [N x )]+C G [Lu 3 -t] 

= f dx [N(x)C(x) + N x (x)D{x) + (w 3 (x) • t) U(x)] , 



(2.44) 



where N, N x and w 3 • t are Lagrange multipliers. The canonical variables evolve with respect to 
the parameter time t via the Hamilton's equations: 



A x = {A x ,n} = 2G 1 



sn 



2G-fiN 



dC 



{ dE x 
A p = {A p ,U} = 2G 1 



N 



dC 



+ N 



5U 



dC 
8E X " 
dC 



2G 7 jiV 
SEP \ dEP 



N 



E x 

E p 

V 
pv 



{E x ,7i} = -2G 7 

{EP,H} = -2G 7 
{ V ,U} = 2G 7 
{ ?? ,^} = -2G 7 



sn 

SA X 

sn 



BG 

-2G 7 { N-^=- + N x 



N 3 



dC 
dEP' 
dD 



dD 
dE*> 



dEP J 



SA, 



- 2G ^\ N dA 




sn 
spv ' 
sn 

Srj 



2G 7 {N x ig + (co 3 .t) 9U] 



-2G 7 



8Pv 



-N- 



K = N x r] + (uj 6 ■ t) 



dP'n J 

= (N x p^y . 



(2.45a) 
(2.45b) 
(2.45c) 

(2.45d) 
(2.45e) 
(2AM) 



The three constraints U = 0, D = and C = together with the above Hamilton's equations are 
completely equivalent to Einstein's equations (with the spherical symmetry imposed). 

It should be noted that the Gauss constraint U = simply yields P v = 0, which is the solu- 
tion to (2.45f). In terms of the canonical variables (A x , A p ,rj; E x , E p , P v ), the Gauss constraint 
(2.42a) involves only the pair (ry; P 71 ), the evolution of which is completely decoupled from that of 



(A x ,A p ;E x ,EP), as the variables (t?;P' ? ) appear only in (2.45e) and (2.45f), while (A x ,A p ; E x , E p ) 



appear only in ( 2.45a )-( 2.45d[ ). In Sec. |II D we will show that the geometry of spacetime is com- 



pletely specified by (A x , A p ; E x , E p ), while (77; P 11 ) correspond solely to the internal degrees of U(l) 



13 Instead of routine calculation, (2.43c I, (2.43d \ and (2.43el can be obtained immediately by knowing the density 
characters of U, D and C (see Appendix |A[) . For example, as D is of density weight 2, according to Footnote |9j 
we have 

{D(x),Cme[N x \} = 2d x N x D + N x d x D, 

which follows 

{C DiB lM x ],Cme[N x ]} = J dxM x {x) \2N X ' (x)D(x) + N x (x)D'(x)} = J dx(M x N x ' - N x M x ')D{x) 
= C DiB [M x N x ' -N X M X '} 
by integration by parts. For |2.43fl, see Appendix [c| for the detailed calculation. 
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gauge. The canonical structure in terms of (A x , A p , 77; E x , E p , P 11 ) thus decouples the U(l) degrees 
of freedom from the geometric/metric ones. 



D. Relation between spacetime metric and Ashtekar variables 



In order to read out the spacetime geometry from the fundamental canonical variables, we need 
a dictionary which translates between the Ashtekar variables and the spacetime metric. 

Take the spherically symmetric 3-dimensional Riemannian space (£, q) and adapting the coordi- 
nates x a to the spherical symmetry: x a = (x, 9, cp), the line element on S is completely characterized 
by two functions A(t,x) and R(t,x): 

da 2 (t) = q ab dx a dx b = A 2 (t, x)dx 2 + R 2 (t, x)dn 2 , (2.46) 

where d£l 2 = d9 2 + sin 2 9d<ft 2 and R is the curvature radius. By the relation qq ab = EfE b 5^ and 



(2.25), we then have 



IE 1 



R 2 



E p = RA. 



(2.47) 



On the other hand, to characterize the line element on the 4-dimensional spacetime (M,g), we 
also need the lapse function N and the shift vector iV a : 



ds' 



g iiv dx p dx v = -N 2 dt 2 + q ab {dx a + N a dt){dx b + N b dt). 



(2.48) 



Because of the spherical symmetry, only the radial component N x of the shift vector survives, and 
both N(t,x) and N x (t,x) depends only on t and x. The extrinsic curvature expressed in terms of 
the lapse and shift functions is given by 



K. 



ab 



1 (8N a dN h dg ab \ _ 1 , , 

+ ~ ~ —. 21 cab l\ ] = (iv a | 6 + iv 6 | a - g ab ) . 



2N \ dx b dx a 



dt 



2N 



For the line element given by (2.46), K ab reads as (also see [23] ) 
E-xx '- 

K ee = - A" " ! R f R - Ft N 



-N~ l A (k - (AA^)') , 
-N^R (l 



Kaa = sin 9K ee , 



and all off-diagonal components vanish. By (2.26), (2.32) and the relation K % a 
—K ab Ej5^ (where e b is the inverse of el), 14 we then have 



(2.49) 

(2.50a) 
(2.50b) 



-K b p i 



A x + (a + P)' = jK* = -jK xx e x 3 
-jK^A- 1 = 7 A" 1 ( A - (AN X )' 



~lK xx 



yP 3 

EP 



A p = 2A p cosa = 2jK^' = -2jK ee e e v = -2jK, ' 



^jKgeR' 1 = 2-fN- 1 (R - R'N 



(2.51a) 



(2.51b) 



It is usually given as K l a = K a b e\ in the literature without the extra negative sign. The sign is merely a convention, 
which flips when the orientation of the parameter time t is reversed. In this paper, we add the extra minus sign 
to be in accord with the Kruskal coordinates. 
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Equations (2.47) and (2.51) make up a dictionary between the metric variables (A, R; N, N x ) 
and the Ashtekar variables (A x , A p ; E x , E p ). More precisely, while (E X ,E P ) gives the intrinsic 
geometry of the leaf St, (A x , A„) tells its extrinsic geometry (i.e. how St is imbedded in M). On 
the other hand, the remaining canonical variables rj and P v correspond solely to the internal U(l) 
degrees of freedom and are completely decoupled from the metric degrees. 



E. Remarks on falloff conditions and boundary terms 

In the full theory, for simplicity, it is assumed that the spatial surface S is a compact 3- 
manifold without boundary. Modifications are required to incorporate boundary terms when E 
is non-compact or has boundary (see [1] for comments and references). In spherically symmetric 
cases, the 1-dimensional manifold I is usually taken to be non-compact or with boundary and 
therefore we have to pay attention to the behavior of the canonical variables at boundary (or at 
infinity if / is non-compact). For the maximally extended Schwarzschild spacetime, which has two 
(left and right) infinities, the falloff conditions of the canonical variables and the corresponding 
boundary terms have been studied in [23J. 

Let (x, y, z) be a system of coordinates on S which is asymptotically Cartesian. Such a system 
is related to a spherical system of coordinates (r, 9, <p) through the standard flat space formulae, 
i.e. (x,y,z) = (r sin cos 0, rsin^sin^, rcosQ). Consistency then demands that, at r — > ±oo, the 
canonical variables have to satisfy the falloff conditions: 

A(t,r) = l + 2GM±(t)\r\- 1 + 0°°(\r\- {1+e) ), (2.52a) 

R{t,r) = \r\ +0°°(\r\- e ), (2.52b) 

P A (t,r) = 0°°(\rD, (2.52c) 

P R (t,r) = 0°°(|r|-( 1+£ )), (2.52d) 

where Pa and Pr are conjugate momenta of A and R, respectively, and f(x a ) = 00 (r _ri ) means 
that / falls off like r~ n , f a like r ~^ l \ and so on for higher spatial derivatives. Meanwhile, the 
Lagrange multipliers have to satisfy 

N(t,r) = N±(t) + 0°°{\r\- e ), (2.53a) 
N x (t,r) = 0°°{\r\- e ). (2.53b) 

It turns out M± are the ADM energy measured at right and left infinities (i.e. the Schwarzschild 
mass) as the boundary terms, and N± correspond to the Lagrange multipliers on the boundary. 
Our main purpose is to study the Kruskal extension of the Schwarzschild black hole (and the 



loop quantum corrections on it). In Sec. Ill, we will see that the Schwarzschild solution in terms 
of Kruskal coordinates leads to N± =0. In this particular situation, the boundary terms vanish 
and no furthermore modifications have to be taken into account. 



III. MAXIMALLY EXTENDED SCHWARZSCHILD SPACETIME 

In order to study the interior and exterior of the Schwarzschild black hole at the same time in 
Ashtekar 's formalism, we need a foliation and a system of coordinates for St, such that St includes 
both interior and exterior and all dynamical variables are well-behaved everywhere on Sj (except 
the singularity). The standard Schwarzschild solution in terms of Schwarzschild coordinates is 
unsuitable for this purpose, as some of the metric components are ill-behaved across the horizon. 
The appropriate coordinates are the Kruskal coordinates, which have the advantage that they 
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cover the entire spacetime manifold of the maximally extended Schwarzschild solution and are 
well-behaved everywhere except the singularity |26t [2"7] . The maximally extended Schwarzschild 
spacetime is also called the Kruskal extension of the Schwarzschild solution. 



A. Kruskal extension of the Schwarzschild black hole 

The Schwarzschild metric is given by 

, „ / 2GM\ 9 ( 2GM\~ X 9 9 ~ , . 

ds 2 = -(l _jdr 2 +fl —J dR 2 + R 2 dtt 2 , (3.1) 

where T is the Killing time and R is the curvature radius. Kruskal coordinates are defined, from 
the Schwarzschild coordinates (T, R, 6, (ft), by transforming (T, R) to the new coordinates (t, r) via: 



R 



2GM 

for the exterior region R > 2GM, and 

1/2 



^"^(Jm)' < 3 ' 2 *» 
1 )" 2e ' i/4GMcosh (^)' (3 ' 2b) 



for the interior region < R < 2GM. In the Kruskal coordinates (t,x,6,4>), the metric of the 
extended Schwarzschild spacetime is given by 

dg2 = 32GHP e _ R/2GM 2 + ^ + R2dn2 ^ (3 4) 

R 

where the curvature radius R = R(t, x) is a function of t and x via 



or more explicitly 

R{t,x) = 2GM 
with W(x) being the Lambert W function. 15 



1 + W 



x 2 - t 2 



(3.6) 



The Lambert W function, also called the Omega function or product logarithm, is the inverse function of f(W) 
We w . The derivative of W is given by 

W(x) 



W'{x) 



x[l + W{x)\ 



17 





FIG. 1: The Kruskal diagram — the Schwarzschild solution in Kruskal coordinates x and t (with 9 and <fi 
suppressed). Contours of constant R are hyperbolae as shown (smaller R in lighter color/shade); particularly, 
t 2 = x 2 corresponds to the event horizons (R = 2GM) and t 2 — x 2 — 1 corresponds to the black and white 
hole singularities (R — 0). Contours of constant T are straight lines passing the origin (not shown except 
for the event horizons, of which t = —x corresponds to T — — oo and t = x to T — +oo). Four regions are 
also indicated. The Kruskal diagram can be conformally deformed into the Penrose diagram as in Fig. [3] (a) 
and (c). 



The transformation between Schwarzschild coordinates and Kruskal coordinates is defined for 
R > 0, R 7^ 2GM, and — oo < T < oo, which is the range for which the Schwarzschild coordi- 
nates make sense. However, the Kruskal coordinates can be extended beyond the range of the 
Schwarzschild coordinates and the allowed values are — oo < x < oo and t 2 — x 2 < 1. The maxi- 
mally extended Schwarzschild spacetime can be divided into four regions as depicted in the Kruskal 
diagram (with 6 and (f> suppressed) in Fig. [TJ The four regions are separated by event horizons 
at t 2 — x 2 = (i.e. R = 2GM). Region I is the exterior region described by (3.2); Region II is 
the black hole interior described by (3.3); Region III is the parallel exterior region described by 
(3.2) with (t,x) replaced by (— t, —x); and Region IV is the white hole interior described by (3.3) 
with (t,x) replaced by (— t, —x). The extended spacetime has two singularities at t 2 — x 2 = 1 (i.e. 
R = 0) for the black and white holes respectively. 



B. Classical canonical solution 



To have a canonical solution corresponding to the Kruskal extension of the Schwarzschild black 
hole, we adapt the foliation in accordance with the Kruskal coordinates: that is, let the leaves St 
be the horizontal lines in the Kruskal diagram and the vertical coordinate t be the parameter time 
of the foliation. 



Comparing (3.4) with (2.48), we then have 



N(t,x) 
N x {t,x) 



R 



0; 



(3.7a) 
(3.7b) 
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comparing with (2.46), we have 



A(t,x) 



32G 3 M 3 
R 



-R/AGM 



(3.8) 



and R given in (3.4) is identical to the curvature radius R defined in (2.46). Equations (2.47) and 



(2.51) then yield 



E x (t,x) 
E p (t,x) 

A x (t,x) 
AJt,x) 



R 2 , 



RA = V / 32G 3 M 3 R 1 / 2 e ~R/^GM 



jN^A 



1 1 
-7 I — + 



2R AGM 



R, 



2-yN~ 1 R = , 27 R l ' 2 R e R / AGM . 
V32GHP 



(3.9a) 
(3.9b) 

(3.9c) 
(3.9d) 



It is tedious but routine to show that (3.9) with (3.6) satisfies D = and C = given by 



(2.42b) and (2.42c) and the Hamilton's equations (2.45) with the lapse and shift functions given 
by (3. 7). 16 Thus, (3.9) is the solution in terms of Ashtekar variables for the Kruskal extension of 



the Schwarzschild spacetime. 



It should be noted that (3.7b|) trivially satisfies (|2.53b|). Meanwhil e, (|3.7 a[) together with (|2.52b|) 

1-172 



leads to N(t, r) = O c 



-|r[> 



which implies that N± (t) 



terms vanish. Furthermore, by A(t , r) 
0°°(|r| 1+e ), which satisfies (|2.52a|) with M±(t) 



in (2.53a) and thus the boundary 
A(t, x)\dx/dr\, we can compute A(t, r) = 1 + 2GM±\r\~ 1 + 



M. 



IV. LOOP QUANTIZATION OF SPHERICALLY SYMMETRIC MIDISUPERSPACES 



As the Ashtekar's formalism for the spherically symmetric spacetime has been formulated for 
classical theory and the classical canonical solution has been obtained for the maximally extended 
Schwarzschild spacetime in the previous sections, we now turn our attention to loop quantum 
theory of the spherically symmetric midisuperspace. The framework has been developed in |18H20| . 
Kinematics of the quantum theory is well understood, but much remains to be done for dynamics. 
Our intention here is not to construct a rigorous loop quantum theory but rather to review and 



bring out conceptual and technical issues. In Sec. IV A, we recapitulate the ideas of previous works 



for kinematics but meanwhile suggest modifying some details. In Sec. IV B we discuss the dynamics 
and propose a new strategy to construct a graph-preserving Hamiltonian constraint operator. In 



Sec. IV C we propound a new quantization scheme inspired by the "improved dynamics" in LQC 
and show that it is a more sophisticated and sensible way to impose the fundamental discreteness 
of LQG when the remnant diffeomorphism is concerned. 



A. Loop representation 

Quantization in the loop representation is based on cylindrical functions of connections through 
holonomies. Let g be a graph on the spatial manifold I composed of a set of edges {e} and points 
{p}. To begin with, we keep the graph g generic: edges of g may or may not overlap with one 
another; points are not necessarily the endpoints of edges and may or may not intersect edges. 



This has been explicitly verified by the authors with the help of Mathematica. 
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When a point intersects an edge, the intersection point is called a vertex. 17 The vector space of 
cylindrical functions with support on a given graph g is denoted as CyL, an element of which, in 
terms of the connection variables A x , A p and rj, is given by 

%(A) ee (A\%) 

= ^2 II ex P C~Y J A x (x)dxj J| exp ( ^ApipU exp (in p r/(p)) , (4.1) 

which is an almost periodic function of ^ f A x {x)dx and ^A p (p) as well as a periodic function of rj. 
In summation, edge labels k e £ M and point labels /i„ G M and n v 6 Z all run over a finite number 
of values. Note that Ar is a scalar density of weight 1 and is integrated along an edge to yield 
(radial) holonomy, while A p and rj are scalars of weight and the cylindrical function is via their 
"point holonomies". As far as cylindrical functions are concerned, an arbitrary graph can always 
be regarded as a set of non-overlapping edges and the union of all endpoints of the edges (thus, all 
points are now vertices and henceforth we rename \i v and Tip CIS f-ly and n v ). 18 Allowing edge and 
point labels to be zero, any cylindrical function can be viewed with support on an (irregular and 
finite) 1-dimensional lattice in /. 

Remark. The vertex labels are integers as n takes values in a circle S 1 (i.e. n = r] + 2ir). The 
other two labels k e and jj, v , on the other hand, are real numbers. Contrast to the full theory, in 
which the edge labels of spin networks are quantized to half integers, the edge labels k e here are 
not quantized. This is because the SU(2) internal gauge group is reduced to £7(1) in the spherical 
symmetry-reduced theory and the weight of irreducible representations of U(l) is continuous. In 
previous works [18H20] . the holonomy := exp \i f e A x (x)dx is regarded as an 5C/(2)-valued T3- 
holonomy and accordingly k e are quantized to integers. In this paper, we adopt a different approach 
by directly considering cylindrical functions of the symmetry-reduced connection variables from the 
outset without adhering to any SU(2) structure inherited from the full theory. Our approach is 
essentially in the same spirit of that in LQC and its merit will be clear later as it makes possible 



the graph-preserving Hamiltonian constraint operator in Sec. IV B which in turn facilitates the 



improved quantization scheme in Sec. IV C 



The vector space of all cylindrical functions is denoted as Cyl, which is the projective limit of 
Cyl ff for all g. The Cauchy completion of Cyl is the kinematical Hilbert space H, which is spanned 
by an orthonormal basis of (symmetry-reduced) spin network states \g; k e , fi v ,n v ): 

(A\g;k e ,ji v ,n v ) = exp f ^ / ' A x (x)dx\ exp ( ^ A p (v) J exp (in v rj{v)) , (4.2) 

where v S v(g) are vertices of g and labels {k e , ji v , n v } are called the coloring of the spin network. 
A spin network defined on g can be regarded as a spin network with support on a larger graph 
g D g by assigning trivial labels to the edges and vertices which are not in g. For any two graphs 
g and g' , let g = gL) g' and we have 

(g; k e , fi v ,n v \g ; k e , fi v ,n v ) = (g; k e , ji v ,n v \g; k e , jx v ,n v ) 

= n h e ,K II 6 ^yjn v ,n' v - (4.3) 

eSe(§) vev(jj) 



From the point of view of the full theory, the edges e are edges along the radial (d x ) direction, while the points p 
are (collapsed) edges along the homogeneous (de,d^,) directions. The vertex in the full theory corresponds to the 
intersection of a point with edges. 

The overlapped segment is regarded as a new edge whose k e is the sum of edge labels of overlapped edges. If an 
endpoint is not a given point in the first place, it is then considered as a point with /i p = n p = 0; on the other 
hand, a dangling point is considered as an endpoint (a conjunction point) of two adjacent "zero edges" with k e = 0. 
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It should be noted that even though k e and fi v take continuous values, S^k' and S^^i are Kro- 
necker deltas rather than Dirac delta functions, as (4.1) takes the form of summation instead of 



integral. The kinematical Hilbert space H is spanned by Tig, the subspace of H with support on 
g, and each U a is given by L 2 (R l £ o hr , S e \ k eBohr ) <gi L 2 {R l ^ ohr , dM/^Bohr) ® where M B ohr is the 
Bohr compactification of R and |e| and |u| are numbers of edges and vertices of g. 

The action of the flux operators that quantize the momenta E x , E p and P v are given by (with 
the Planck length £p\ := V Gh ) 



E x (x)\g;k e ,fj, v ,n v ) = 7^ e+(a) ^~ e ^ \g;k e , fi v ,n v ), 



(fx-P^a;)^; k e ,fx v ,n v ) 



v&v(g)nl 



^e? 
■uS'u(g)nX 



(4.4a) 
(4.4b) 

(4.4c) 



where e^{x) are two edges or two parts of a single edge meeting at x and I C / is an arbitrary 
region of I. 19 Note that E p and are of density weight 1 and thus their corresponding operators 
are to be smeared over a region of /, whereas E x is of density weight and its corresponding 
operator is not smeared. 

Consequently, after regularization procedure as in [19], the volume operator corresponding to 
the volume V(T) = Att j x dxy/\E x \ E p of a region I C / is given by 20 



V{X) = Air J dx\E p \\l\E x \ , 



where E p (x) is the distribution-valued operator: 



E p (x)\g;k e ,fJ, v ,n v ) = j£p l ^ S(x - v) fj, v \g;k e , fj, v ,n v ). 



(4.5) 



(4.6) 



The volume operator then has \g; k e , fi v , n v ) as its eigenstates: 



V(l)\g;k e ,Hv,n v ) = 4vr 7 3/2 4i ^ 1^1 

v&v(g)rd 



\g, k ej fi v , n v ) 



(4.7) 



The Gauss constraint with (2.20) is classically given by 



C G [\} = ^J i dx\(E x ' + pv)^0, 



and quantized to 



f 2 
2G 



E ^) {k e +(v) ~ K-(v) + 2n„) \g; k e ,fiv,n v ) = 0, 



(4.8) 



(4.9) 



vGv(g) 



19 If x does not intersect any edges, we have k e +r x \ = k e -r x \ =0. If a; is the right (resp. left) endpoint of e but does 
not touch another edge, we have k e +^ = (resp. k e -r x \ = 0) 



Note that both A p (as well as A p ) and E p are defined to be non-negative as in ( 2. 16 1 , and thus only /j, v > is allowed. 
However, it is technically easier to allow all values \i v £ H and in the end require physical states to be symmetric 
under the large gauge transformation fi v i-» —fi v (see [9 ). Classically, flipping the signs of A p and E p simultaneously 
yields the equivalent solution, as the constraints (2.421 admits the large gauge (A P ,E P ) — > (— A p ,— E p ). 
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which is solved exactly by 



n. 



e+(v) ke~ (v) 



)) 



(4.10) 



Substituting (4.10) back to (4.11), we obtain the general form of the U(l) gauge invariant spin 
networks \g;k e ,fi v ) £ Uf'. 



(A\g;k e ,fi v 



e£e(g) 

n 



cxp 



eSe(g) 



ik e 
~2~ 

2/fcjp 



L (A ' 



+ rf)dx ) ] [ exp 

vev(g) 

ifj, 



AJv 



A x dx ) ] [ exp 

v£v(g) 



(4.11a) 
(4.11b) 



which are cylindrical functions of the U(l) invariant connections A x and A p . In accordance with 



the canonical relations (2.41), the flux operators in (4.4) acting on the f7(l)-invariant spin networks 
are given as 



E x (x)\g;k e ,fi v ) 



+ h 



e (x) 



9'i kei t^v) > 



J dxE p (x)\g;k e ,fj, v ) = 7^ ^ fj, v \g;k e , /j, v ), 



v£v(g)nX 



dxP v (x)\g;k e ,fi v ) = 0, 



(4.12a) 
(4.12b) 

(4.12c) 



where the last equation is equivalent to the Gauss constraint. Similarly, the volume operator acting 
on the [/(l)-invariant spin networks is given by 



V(l)\g;k e ,n v ) = 4ir 1 3 / 2 e 3 Pi £ 

vEv(g)riZ 



\g] k e , fi v ) . 



(4.13) 



The next is to impose the diffeomorphism constraint Coif^A^] ~ 0. Following the strategy used 
in the full theory (see Section 6.2 of pQ for the details), the procedure of group averaging can be 
applied in the similar fashion to obtain diffeomorphism invariant states. Although one begins with 
a state in Cyl, after group averaging, the resulting diffeomorphism invariant state is distributional 
and belongs to Cyl*, the algebraic dual of Cyl. The subspace of all elements of Cyl* which are 
invariant under diffeomorphism is denoted as Cyl*j iff . In the end, every element (^| 6 Cyl*j iff can 
be uniquely decomposed as 



(*i=E(%ii> 

Id] 



(4.14) 



where [g] runs through the diffeomorphism classes of graphs. If two graphs g and g' are not in the 
same diffeomorphism class, (^[ g ] \ and are orthogonal to each other, i.e. (\Pf g i = 0. The 

basis states of "Hoiffi the Cauchy completion of Cyl*) iff , are first labeled by the diffeomorphism class 
of graphs and then distinguished only by colorings of edges and vertices. The colorings labeled 
for spin networks in T~L are not necessarily orthonormal in "HDiffj due to the nontrivial action of 
the discrete graph symmetry group. 21 The states of an orthonormal basis of %Diff ar e called 



21 The graph symmetry group in 1-dimensional manifolds is much simpler than that in 3-dimensional manifolds, yet 
non-triviality still remains. 
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(symmetry-reduced) spin-knot states or s-knot states (by adopting the same name used in the full 
theory), each of which is specified by a graph diffeomorphism class [g] and a discrete coloring of 
edges and vertices, which is different from the coloring of spin networks. 

Disregarding the technicalities due to the graph symmetry group, one can view a spin-knot as a 
spin network state by dismissing any reference to its localization in I. Similar to the interpretation 
in the full theory (see Figure 6.8 of [2]), a spin-knot can be abstractly interpreted as an ensemble of 



chunks of volumes, given in (4.13), arranged in a radial order with adjacent surfaces of areas, given 
in (4.12a). 22 A spin-knot state possesses the information about the volumes and the areas that 
separate these volumes, but any information of localization of the chunks and adjacent surfaces 
is irrelevant. However, contrary to the full theory, where the eigenvalues of volumes and areas 
are quantized, the volume and area eigenvalues are only "partially" quantized: that is, they take 
continuous values but the values are Bohr compactified. This is a consequence that we reduce the 
internal SU(2) gauge to U(l) from the outset. 

B. Hamiltonian constraint 

The most complicated and difficult part is to quantize the Hamiltonian constraint. In the full 



theory, a standard procedure is established to quantize (2.13) by taking care of three complications 
[1]: first, curvature components F* b are expressed in terms of holonomies which can be directly pro- 
moted to operators; second, the operator corresponding to inverse volume element e _1 is obtained 
from a commutator between holonomies and the volume operator a la Thiemann's trick; third, 
the Lorentzian part of C[N] (which involves extrinsic curvature K^) is expressed as a commutator 
between the total volume and the Euclidian part of C[iV]. This procedure can also be adapted to 
symmetry-reduced models (see [3HH] for the case of LQC). For the spherically symmetric models, 
the resultant Hamiltonian constraint operator has been derived and given in Equation (41) of [20J, 
which essentially is constructed in terms of the volume operator V{v) and the holonomy operators 
acting on the vertices of spin networks; that is 

C[N]=J2N(v)C v , (4.15) 

V 

and C v is expressed in terms of V(v) and the holonomy operators 

(I Mv)±e \ 
- / A x (x)dxrz , (4.16a) 
^ Jx(v) J 

h e {v,S) = exp(^A p (v)T r ) , (4.16b) 



h<p(v,5) = exp Ap(v)T 2 ^j , (4.16c) 

where e is the coordinate length of the radial edge starting at the vertex v and ending at the 
coordinate x(v) ± e and the sign ± denotes the orientation of the radial edge running to the right 
(+) or left (— ) of the vertex v. The parameters e and 5 are used for regularization. 



The SU(2) values arising from the exponentials of T3,T\',T2> in (4.16) will be traced out in 



the end and thus yield only combinatorial numbers. Apart from the combinatorial numbers, the 

22 More precisely, each of the chunks of volumes takes the shape of a spherical shell, as the adjacent surface takes 
the shape of a sphere. 
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operator C v essentially involves only the volume operator V(v) and the two holonomy operators 
for the radial and point holonomies: 



h 



exp 



K(v) 



x(v) 
f i5 A f > 

exp I — Ap{v t 



A x (x)dx 



(4.17a) 
(4.17b) 



where we have replace A x with A x for h%±(v) as from now on we will restrict ourselves to the U(l)- 
invariant spin networks. When acting on a vertex v of a [/(l)-invariant spin network \g;k e ,fi v ), 
the aforementioned operators only interfere with the vertex v, its two neighboring vertices and its 
two connecting edges. We therefore drop all other labels and denote 



\n-,k-,fj,,k + ,fi + ) 



k- v k-> 



The actions of the operators are then given explicitly by 



and 



V(v)\fi-,k-,fi,k + ,fj, + ) 
h 5 (v)\^-,k-,fi,k + ,fi + ) 



K t+ (v)\fJ,-,k-, n,k+, /J,+) 



h% (v)\fjt-,k-,n,k+,ii+) 



4vr 7 3/2 4ll/^l/ 



k+ + k. 



\fi-,k-,fi,k + ,fi+), 



\fi-,k-,/j, + S, k+,(M+), 



if e > 

k-,(i,k + +l,fi + ), if e = 

fe++l 3 0, fc+,//4.}, if e < 

if e > 

0, — 1, if e < 



Ax + 
Ai + 

Ax_ 



(4.18) 

(4.19) 
(4.20) 

(4.21a) 
(4.21b) 



where Ax± is the coordinate distance from the vertex v labeled by fi to the neighboring vertex 
labeled by fjL±. When e < |Aa;±|, the endpoint of the radial holonomy does not fall on vertices 
of the original state and thus creates a new vertex in between labeled by 0. Note that V{v) is 
hermitian, while h s Jv) and h%±{v) are unitary. 



The Hamiltonian constraint operator implemented with (4.21) changes the graph of a given 



spin network state by creating new vertices and thus makes the dynamics rather complicated. 
Alternatively, we can construct a graph-preserving Hamiltonian constraint operator, which changes 
only the coloring of spin network states, by devising a different regularization scheme modified from 



(4.17a): 



K,±W) = ex P 



IE 



A x (x)dx , 



(4.22) 



2 |Ax±| 7 e ±(„) 

which, when acting on spin network states, gives 

h% k-, (J,, k+, = k-, /x, k + +e |Ax + | _1 , (4.23a) 

h% (J,, k+, (/,+) = k^ — e |Ax_| _1 , fi, k + , /x + ). (4.23b) 

The regularization variable e is now introduced to change the edge label of the connecting edge, 



instead of attaching a new edge starting from the vertex. The action of (4.23) can be viewed as 
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remark in Sec. IV A). 



evenly stretching the new edge in (4.21) to the whole extent of the connecting edge e =t (u). The 
quantization scheme (4.23) is possible only if we allow edge labels k e to be real numbers (see the 



In the full theory, the Hamiltonian constraint operator, albeit defined upon spin networks with 
reference to coordinates in the first place, turns out to be well-defined and independent of the 
regularization at the level for diffeomorphism invariant states, i.e. spin-knots, as a consequence of 
the intimate interplay between diffeomorphism invariance and short-scale discreteness [U [2] . The 
remarkable feature of LQG that the dependence of the regulating parameter disappears at the 
level of spin-knots is no longer the case in symmetry-reduced theories, because the diffeomorphism 
invariance as well as the SU(2) gauge invariance of the full theory are partially (in the case of 
spherically symmetric models) or completely (in the case of LQC) broken, and as a result the 
short-scale geometry is only "partially" discrete (in the sense of Bohr compactification) . That 



said, (4.20) and either of (4.21) and (4.23) are well-defined only upon spin networks, and cannot be 



carried over for spin-knots. On the other hand, the Hamiltonian constraint operator is ill-defined 
in the limit e, <5 — > 0, because of the very nature of loop representation that connection operators 
do not exist. (However, see Section 5.3 of [20] for more comments on and alternative approaches 
for regularization issues and anomalies.) 

In order to faithfully manifest the discreteness of quantum geometry of LQG in symmetry- 
reduced models, we follow the strategy devised in LQC: keep the regulating parameters e and S 
finite by hand to imprint the fundamental discreteness of LQG in a sophisticated manner such that 
the Hamiltonian constraint is diffeomorphism-invariant. The next subsection is devoted to issues 
of this strategy. 

C. Improved dynamics 



In LQC, to impose the fundamental discreteness, if the regulating parameters are simply pre- 
scribed to be small constant, it has been demonstrated that this prescription (i.e. /icrscheme) leads 
to wrong semiclassical behavior, due to the underlying fact that the prescription is not independent 
of the choice of the comoving finite sized cell to which the spatial integration is restricted to make the 
Hamiltonian finite. To fix this problem, a more sophisticated quantization scheme of "improved" 
dynamics (i.e. /x-scheme or, more precisely, //'-scheme for anisotropic models) was formulated, in 
which the regulating parameters are taken to be adaptive discreteness variables. The quantization 
scheme of improved dynamics was first developed in [T3] and later generalized and elaborated on 
for many other models (see |15|ll6j for Bianchi I models and |13|ll7j for Kantowski-Sachs models). 

In our case of spherically symmetric midisuperspaces, the same problem of wrong semiclassical 
behavior is anticipated if we simply take e and 5 to be constant, because the 1-dimensional dif- 



feomorphism invariance will be broken as both (4.21) and (4.23) depend on the radial coordinate 
(either implicitly or explicitly). Therefore, we have to formulate a new quantization scheme which 
resolves coordinate dependence. 



One might argue that the quantization scheme (4.231 may not be legitimate, as the connecting edge of a given 



spin network could be quite long and the holonomy along the long edge no longer well approximates the curvature 
strength at the vertex, thus, leading to breakdown of approaching classical dynamics in the classical regime. 
This, however, is not a problem. A spin network state which represents a smooth classical geometry is called a 
weave state [2], in which the physical lengths of edges are supposed to be so small that the granular structure of 
space is completely negligible as far as attainable technology is concerned. The quantization scheme ( |4.23[ ) indeed 
approximates the classical dynamics very well when acting on the weave state. 

If one formally take the limit e, S — > on the Hamiltonian constraint operator, it will give rise to the Wheeler- 
DeWitt quantization, which is essentially different from the loop quantization (see |3T] for the case of LQC). 
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Ashtekar's formalism for spherically symmetric models has a direct correspondence to that for 
the Kantowski-Sachs spacetime if the homogeneity is formally imposed. It is not surprising that the 
improved quantization scheme for spherically symmetric models bears close resemblance to that 
for the Kantowski-Sachs model as formulated in [17]. Following the same ideas (of the //'-scheme) 
in Appendix B of [T7j, areas of rectilinear loops of holonomy have to be shrunk to the area gap 
obtained in the full theory of LQG. Let p, x and jl p be the coordinate lengths of the holonomy paths 
in the radial and homogeneous directions, it is required (by the so-called /i'-scheme) that 



tixUpAR = 2A, 
JxlR 2 = A, 



(4.24a) 
(4.24b) 



where A = is the area gap of LQG and £ = 2v3 ir for the standard choice (but other choices 

are also possible). By (2.47), the discreteness variables are solved as 



t'.i 



\E X \A 



= 2 



A 



A 

\E X \ 



R 



(4.25) 



Note that we have shrunk the area in (4.24a) to 2A while the area in (4.24b) to A. The difference 



of the factor 2 gives rise to the extra factor 2 for ji x in (4.25). This factor is demanded by 



the consistency of diffeomorphism and Hamiltonian constraints, as we will see in Sec. VA This 
extra factor was not introduced in [T7j for the Kantowski-Sachs models, since the diffeomorphism 
invariance is broken (and reduced to scaling invariance) by homogeneity and thus the ratio between 
the areas of (|4.24ah and (|4.24bl) is not fixed. 25 



Replacing 5 in (4.20) and e in (4.23) with ±[J, P and ±[i x prescribed in (4.25), respectively, we 
can accordingly construct the new holonomies h^ p (v) and 
network states, are given by 



v), which, when acting on spin 



h^(v)\fi-,k-,fi,k+,n+) 




(4.26a) 
(4.26b) 
(4.26c) 



in accordance with the actions of flux operators given by (4.12 ). 26 Note that the new holonomies h p p 



and h^ x ± remain unitary, but (/i^_j_) _1 / h~± while (h p p )~ l = h p ^ p . The improved quantization 
scheme with (4.26) is now independent of coordinates. However, a minor complication arises from 
the improved quantization scheme: the regulating factor 5~ 2 appearing in the leading factor of 
Equation (41) in [20j is now replaced by p,~ 2 , which is given by \E X \/A. As a result, the flux operator 



Heuristically, the factor 2 in (4.24a I can be understood as a tradeoff that (4.24a I actually represents two holonomy 
loops in (d x ,de) and (dxyd^) directions, but the two degrees of 9 and <f> are collapsed to one by the spherical 
symmetry. 

The holonomy operators implemented in (4.261 are reminiscent of Equation (3.15) in [TS] for LQC of Bianchi 
I models. The techniques of algebraic simplification in [16] by introducing new variables could be carried over 
to handle the complicated operators in a manageable fashion for further rigorous developments of the quantum 
theory. By the way, the same reasoning following Equation (3.19) of [IB] can be used to argue that the operators 
in (4.26 1 are well-defined even though the denominators of the radicands could be zero. 
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E x {v) acting at the vertex v is also needed to express this factor. Once this minor issue is taken 
care of, the resultant Hamiltonian constraint operator is cast in terms of h^ p (v), h x ^(v), V(v) 
and E x {v). Therefore, information of localization of spin-network states is completely irrelevant 
and the dynamics dictated by the Hamiltonian constraint operator is diffeomorphism invariant. 

Although the improved quantization scheme respects the diffeomorphism invariance, it gives 
rise to the problematic feature: [C v , C v >] ^ if v and v' are adjacent to each other. Consequently, 
given different smearing functions, the Hamiltonian constraint operators do not commute, i.e. 
[C[iV], C[M]] 7^ 0, even restricted to the sector of diffeomorphism invariant graphs. The physical 
implication is that the resulting dynamics depends on foliation of the spacetime. This problem 
results from the imposition of the fundamental discreteness, which breaks down the intimate match- 
ing between the diffeomorphism and Hamiltonian constraints in the full theory. However, inspired 
by the ideas for in |24|I25| for LQC, one might speculate that a linear superposition of the Hamilto- 
nian constraint operators in generic j representations might conspire to give the new Hamiltonian 
constraint operator by which [C v ,C v r] = 0. That is, if we consider all generic j representations 
rather than a fixed one (particularly, j = 1/2), the notion of permissibility of the classical regulator 
could be restored (at least at large scale). As studied in [25J, inclusion of higher j representations 
corresponds to higher order holonomy corrections. While rigorous construction of the quantum 
theory is rather difficult, without further ado, in Sec. |VJ we will investigate the ramifications of 
the improved dynamics with higher order holonomy corrections at the level of heuristic effective 
dynamics. 



V. HEURISTIC EFFECTIVE DYNAMICS 

In order to see the ramifications of the loop quantum corrections without going into the de- 
tailed construction of the quantum theory, one can study the effective dynamics at the heuristic 
level by viewing the dynamics as classical but governed by the new "holonomized" Hamiltonian 



constraint, which, to capture the loop quantum corrections, is modified from (2.42c) by replacing 
the connection variables with the holonomized ones: 

A x -> sm W* A *> = e _ _? ; (51a) 

[i x 1i[ix 

P ftp 2ifi P 

This prescription incorporates "holonomy corrections" but ignores corrections due to Thiemann's 
trick (particularly, the "inverse triad corrections"). It was argued in [31] that, for LQC, the cor- 
rections due to Thiemann's trick is negligible on physical grounds. Various studies on LQC have 
suggested (and verified in some particular models) that the heuristic effective dynamics with holon- 
omized Hamiltonian constraint indeed gives a very good approximation to the quantum evolution 
for the states which are semiclassical at large scale [14|, [32Tl34j . For the spherically symmetric 
models, it is also expected that the heuristic effective dynamics can provide us considerable insight 
about the quantum evolution. 



A. Higher order holonomy corrections 



Following the ideas in 



5], the prescription in (5.1) can be refined in a more elaborate way 



to include higher order holonomy corrections. This to replace the connection variables with the 
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FIG. 2: The higher order holonomized connections A x n p as functions of A x ,p (for given p x ,p)- The periodic 
functions A x n J, (shown with n = 0, 1, 2, 3) asymptote to the triangle wave A x °^ . A x °^} agrees with A XtP for 
-7r/2 < fi x ,pA X)P < 7r/2. 



nth order holonomized ones: 



(2fc)! 



^^22^(^)2(2^ + 1) 



sin 



((i x Ax)) 



2fc+l 



4 n) : =-E 



(2k)\ 



k=0 



2 2k {kl) 2 (2k + l) 



sin 



(5.2a) 
(5.2b) 



which can be made arbitrarily close to A x an A p (as n — > oo) for —tt/2 < fi x ,pA x , P < vr/2 but remain 
functions of the holonomies ht^ x '■= ex.p(±ifi x A x ) and h^^ p := ex.p(±iji p A p ) as well the discrete- 
ness variables fl x and p> p (see Fig. [2]). The higher order corrections correspond to higher powers of 
s'm(fl x A x ) and sin(p, p A p ), which might be understood as a result of generic j representations for 
holonomies in the Hamiltonian constraint operator in the quantum theory [25j . 

It is noteworthy that the density weig ht of A x n) is the same as that of A x , and A p the same as 
A p (see Appendix |a|). Substituting (5.2) into (2.42c), we obtain the Hamiltonian constraint with 
the nth order holonomy corrections: 



(j{n) ._ 



EP 



2G \ ^/\E^\ 4 7 2 VP 

-sgn(^) VI + 



E" . rs JW\A p n >A 
-sgn(^)^ 1 1 ' 



X 



r 



(E 3 



xl\2 



(E 



p\2 



ly/\E*\EP 



+ sgn(£ a 



\fW\ EX 

EP 



(5.3) 



which has the sa me den sity weight as the Hamiltonian constraint C . Consequently, the constraint 
algebra given by (2.43e) remains unchanged if C[N] is replaced by C^[N] := J f dxN(x)C^(x): 



C^[N],C DiS [N x ]} = -C (n) [N x N'} = -C^ n) [C N ,N], 



(5.4) 



indicating that the diffeomorphism invariance is still respected regardless of the order n of holonomy 
corrections. 



28 



Unfortunately, as the imposition of fundamental discreteness breaks down the intimate matching 
between diffeomorphism and Hamiltonian constraints, the constraint algebra given by (2.43f) does 
not hold any more KC[N] is replaced by C^[N]. This can be seen as the Poisson brackets between 
the triad variables and the holonomized connections are no longer the same as before in (2.41 ) but 
instead given by: 



{A p n \x),E»(y)} 
{A x m \x),A p n \y)} 



2Gj5(x - y) cos(fx x A x ) ^ 2 ^(ki) 2 ( sin (^^)) 2fc , 



k=0 
it 



2G-/5(x - y) cos(n p A p ) ^ 2 ^(k\) 2 (^(Mp^)) 2 * 



(5.5a) 
(5.5b) 



2G 1 5(x -y)J2H WTfm (^(fix**)) 

~cos(fi x A x ) ( - - sm(fl p A p )\ dftp 



A x cos(n x A x ) 



2fe / . ,_ j s\2l 

(sm(fipAp)) 



{AW(x),E>>(y)} 
{A x m \x),A x n \y)} 



cos(flpAp) 

fix 

A p m \x),A p n \y)) = {E*(x),EP(y)} = 0, 



(2l + l)fipJ dE* 
sm(fl x A x ) \ dfi x 



{2k + l)fi x J dEP 



(5.5c) 

(5.5d) 
(5.5e) 



where we keep the discreteness variables \x x and \i p as generic functions of E x and E p . 
In the limit n — > oo, the holonomized connections (5.2) take the form of 



^4(0°) 



4(0°) 
p 



7T 

2J1 X 
2ft p 



f h (li x A x ), 



Ja(^ p A 



(5.6a) 
(5.6b) 



where f/\(x) is a triangle wave function of x with period 2tt and /a(0) = 0, /a(±vt/2) = ±1 (see 
Fig. 2). Therefore, A x °°^ and Ap°°^ agree exactly with A x and A p respectively, before the evolution 
approaches the critical conditions fl x A x = ±tt/2 and ft p A p = ±7r/2, which signal occurrence of the 
quantum bounce. Furthermore, by the identity 



^ (2AQ! . 2k 



k=0 



2 2k {k\f 



COS X 



the Poisson brackets in (5.5) become 

{A x °°\x),E x (y)} = 2G 1 5(x-y)f n (fl x A x ) 1 
{A p °°\x),EP(y)} = 2G 1 S(x-y)HflpAp), 

{A x °°\x),A p °°\y)} = 2G 1 5(x-y) 



HfiM [f n {u P Ap)A p -A^ 

-fn(flpAp) (Un x A x )A x -A x ^ 



1 djip 
Ji P dE* 
1 dfi x 
fix dEP 



(5.7) 

(5.8a) 
(5.8b) 

(5.8c) 



and all other brackets vanish, where fn(x) '■= sgn(cosx) = ±1 is a square wave function with 
period 2ir. 
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Before and after the bounce, by (5.6), we have 



p 



and 



±1- 

Pp 



before bounce, 
A x , after bounce, 

before bounce, 
An, after bounce, 



fn(HxA x ) = fn(/J.pA p ) 



+1, before bounce, 
— 1, after bounce, 



(5.9a) 
(5.9b) 

(5.10) 



as fni^xAx) and f n (fi p A p ) flip signs when the evolution jumps over the bounce. 27 Consequently, 
before the bounce, the Poisson brackets in (5.8) are the same as the classical ones in (2.41). On 

(5.11a) 
(5.11b) 

1 dU P 1 dfl x , m 



the other hand, after the bounce, ( |5.8[ ) reads as 

{Ag°)(x),E*(y)} = -2G 7 5(x-y), 
{A p °°Xx),E<>(y)} = -2Gj8(x-y), 

{A x ^(x),A^\y)} = ±2nG 1 8(x-y) 



Hi dE x & BE. 



where we have used the particular prescription (4.25) to yield zero on the right hand side of (5.11c) 
and note that the extra factor 2 for p, x in (4.25) is essential. Therefore, after the bounce, the 
Poisson brackets in (5.11 ) are the same as the classical ones in (2.41 ) if A x and A p are replaced by 



-Ax^ and 



Ap°°\ The extra minus sign indicates that the evolution after the bounce behaves 



as if time was rey ersed. It is remarkable that, before and after the bounce, the diffeomorphism 
constraint (2.42b) can be cast in terms of A x °°^ and A^" 1 as 
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2G 7 



(A' p E p + v'P v - A X E X ') 



' ^(a^'ep + v'Pv-Ax^e* 1 
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2G7 
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2G7 



P X 



4(00) 



E xl 



before bounce, 
after bounce, 



A^) EP + rfP" 



-A 



(00) 



E d 



before bounce, 
after bounce, 



(5.12a) 



(5.12b) 



where the prescription (4.25) is used and again the extra factor 2 is essential. Once more, after 



the bounce, — Ax^ and —A p ^' play the roles of A x and A p as diffeomorphism is concerned. 
Equation (5.12) is used to derive (C6), which leads to 



(bo) 



E 2 

(NM' — MN')-j-j 



C g 



F 2 

(NM'-MN')-r^rf 



(5.13) 



As the evolution approaches the putative singularity, A x and A p increase toward infinity until each of fi x A x and 
flpAp reaches the critical value ±tt/2. Once the critical value is reached, l^l I (resp. \A P I) starts to decrease 
and fn(p-x-Ax) (resp. fn(fi P A p )) flips its sign, thus causing the bounce. It is possible, however, that when p, x A x 
approaches ±7r/2, fJ, p A p approaches with an opposite sign. We will discuss this case in Sec. VB but assume 

fl x A x and fipAp are of the same sign in this subsection. 
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where the symbol = is used to indicate that the identity breaks down briefly during the transition 
period of the bounce, as fn(p>xA x ) and fn(ft p A p ) flip signs at close but slightly different epochs 
during the transition. The breakdown leads to the problem that the two constraints C^ 00 ) = 
and D = are violated after the bounce, implying that the spacetime is no longer invariant 
under spatial diffeomorphism and change of spacetime foliation after the bounce. However, if the 
characteristic parameter used to describe the solution in the classical regime is much larger than 
the Planck unit (for the Schwarzschild spacetime, it is the case that the Schwarzschild mass M 
is much bigger than the Planck mass mp\ := y/h/G), the transition period is so brief compared 
to the characteristic time of the whole spacetime (i.e. GM for the Schwarzschild spacetime) that 
violation of C(°°) = and D = is negligible at large scale. Therefore, if we go to the limit n — > oo, 
the intimate matching between diffeomorphism and Hamiltonian constraints is restored at large 
scale. 

Finally, it is trivial to show 



{c G [\],C {n) [N}} = 0. 



(5.14) 



To summarize, with (5.4), (5.13) and (5.14), the constraints C G [X], C B is[N x ] and C (oo) [iV] yield 



the same constraint algebra as in the classical case (2.43) except for the brief breakdown in the 
transition of the bounce. The brief breakdown in the transition period causes the two constraints 
C(°°) = and D = to be slightly violated after the bounce, whereas the Gauss constraints G = 
remains exact at any epoch and everywhere. This is a consequence of the fact that imposition of 
fundamental discreteness spoils the intimate matching between diffeomorphism and Hamiltonian 
constraints. However, the violation of C^°°^ = and D = can be restored at large scale if we 
go to the limit n — > oo and therefore the heuristic effective dynamics dictated by the modified 
Hamiltonian constraint given by (5.3) with n = oo still provides a reliable description for the large- 
scale physics. Details of the resultant dynamics, however, depend on the foliation of spacetime, 
as by (2.51) the quantities p, x A x and fi p A p , which indicate how strong the quantum correction is, 
read as 



f^xA x — 


-27^ 


n±K xx 


= 2 7 A 


/A 


P-pA p = 


-27^ 
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/a 


which still involve N and N x . 28 











A - (AN*)' 

iVA 
R - R'N X 

Wr ' 



(5.15a) 
(5.15b) 



Until now there is no systematic procedure which leads us to the dynamics of the symmetry- 
reduced theory from that of the full theory of LQG, but symmetry reduction in general can be 



Particularly, if we choose N x = 0, we have 



/'. 



X A X = 27 VA 



1 dA 

A dr' 



A p = 2 7 VA 



1 dR 
Rdr' 



where dr = Ndt and r is the proper time. For the models which admit a foliation in which the spatial slices 
are homogeneous (such as [131 117|). both Ji x A x and p. p A p acquire physical meanings: the former is the Hubble 
rate in (dg,d^) direction (times 27\/A) and the latter is the Hubble rate in x direction. This also explains 



why the extra factor 2 is needed (c/. Footnote 25 1. For generic cases without homogeneity, the dependence on 
the foliation can be interpreted as reflection of the semiclassical traits (spreading, squeezing, etc.) of the weave 
states at the level of effective dynamics; i.e. different choices of foliation correspond to different weave states with 
different semiclassical properties. In Sec. |V C| we choose the foliation in accordance with the Kruskal coordinates 
to investigate the quantum Schwarzschild spacetime, since the Kruskal coordinates give a well-behaved foliation 
as discussed in Sec. |III| Other well-behaved foliations are also possible and they could yield qualitatively distinct 



effective solutions (see comments in the end of Sec. VC I 
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understood on the suppositions that in the full theory there is a regime where the nonsymmetric 
degrees of freedom do not affect too much the dynamics of the symmetric degrees and that the 
states of concern happen to be within such regime. In other words, the symmetric degrees of 
freedom can be treated as "heavy" degrees of freedom while the nonsymmetric degrees treated as 
"light" ones, in the sense of the Born-Oppenheimer approximation. The fact that in spherically 
symmetric theories the constraints violate briefly during the transition period suggests that the 
degrees of freedom of the full theory cannot be separated into spherical (heavy) and non-spherical 
(light) ones in the vicinity of the quantum bounce. That is, in the full theory of LQG, a given 
weave state (coherent state of spin networks) which, after coarse-graining, represents a smooth, 
spherically symmetric space at large scale will eventually manifest granularity of spin networks 
when approaching the quantum bounce. After the bounce, the weave state evolves to become 
smooth and spherically symmetric again if the transition period is very brief. Note that in the 
heuristic effective dynamics it is possible that fn{fixA x ) and f n (fl p A p ) do not flip signs at close 
epochs (or even one of them never flips signs); this situation suggests a different scenario in the 
full theory: approaching the putative singularity, the smooth weave state descends into granular 
spin networks instead of being bounced back to a smooth weave. 



B. Further phenomenological modification 



In Sec. VA, we assume that p. x A x and p> p A p approach the critical value ±7r/2 with the same 
sign as commented in Footnote 27 However, it is possible that when fi x A x approaches ±7r/2, 
JipAp approaches =Fvr/2 with an opposite sign. The extended Schwarzschild spacetime given by 



(3.9) gives an example. When this is the case, (5.11c) is modified to 

1 dfip , 1 dfi x 



{AtHx),A p °°\y)] = T 2nG 1 6(x-y) 
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(5.16) 



and (5.12a) modified to 
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(5.17) 



We do not have |Ai°°^(x), ^4p°°^(y)| = and (5.12b) any more, and consequently (5.13) is no 



longer valid. 



There is a way to fix this problem by modifying (4.25) with inclusion of an extra factor of 



powers of E x . For the first case that [i x A x and [i p A p approach the critical value ±7r/2 with the 
same sign, (4.25) can be generalized to a more generic prescription: 



Ma 
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\E*\ ' 

\E x \ xl 
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for \E X \ > C 2 , 
for \E X \ < C 2 , 



for I < 



1 



(5.18a) 

(5.18b) 
(5.18c) 
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Note that (5.18) reduces to the original proscription (4.25) for I = 0. For the second case that 



/j, x A x and \ipA p approach the critical value ±7r/2 with the opposite sign, (4.25) has to be modified 
as 



Mi 



for \E X \ > C 2 , 

for \E X \ > C 2 , 
for \E X \ < C 2 , 



21-1 EP 
21 
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\E*\ 



\E X \ 
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(5.19a) 

(5.19b) 
(5.19c) 



Here, we introduce a new constant length scale C satisfying A < C 2 -C G 2 M 2 with M being 
the mass constant of the spacetime, assume the bounce takes place when \E X \ < C 2 , and require 



I < 1/2 for (5.18) and I > 1/2 for (5.19) to ensure the convention /j, x , [i p > 0. Note that, in the new 
prescription of (5.18) and (5.19), the dimensions and density weights of \± x and fl p are the same as 



those of the old ones given by (4.25). Therefore, the constraint algebra (5.4) is unchanged with the 



The good features obtained in Sec. 
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the classical solution (3.9) yields 
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new prescription. Furth ermo re, it is not difficult to show that, with (|5.18|), ( |5.11c [) and (|5.1 2b[) ar e 
still satisfied, and with (5.19), (5.16) leads to 



VA (5.13) in particular, are restored. 



and (5.17) leads to 



5.12b). 



When the new prescription (5.19 ) is used for the Schwarzschild spacetime in Kruskal coordinates, 
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(5.20a) 
(5.20b) 



which follows 
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By (3.5), (3.6) and Footnote 15\ we have 
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(5.22) 



which then leads to 
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fipAp 
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(5.23a) 
(5.23b) 



33 



In the limit t 2 



±1 (i.e. R — > 0) towards the (black/while hole) singularity, 



_^ 4t 7 VAGM r21 _ 3/2 

^ x x t 2_ x 2^ ±1 ^/2(4/-2)£ 2; 
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(5.24a) 
(5.24b) 



As the strength of loop quantum corrections are indicated by \x x A x and n P A p , both in (5.24) are 
supposed to blow up when computed with the classical solution. This restricts the value of I to be 



1/2 < I < 3/4. 



(5.25) 



Moreover, one can further fix Z = 3/4 — e with a small positive number e to have (5.21 ) read as 



(J>xA x 



;i + 4e + 0(e 2 )) 1 + 
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fl + 4e), fori?«2GM, 



(5.26) 



so that \x x A x gets bounced almost at the same time as fJ, p A p and thus the breakdown in (5.13) 
during the transition period is optimally mitigated. 



It should be emphasized that, until now, the modification in (5.18) and (5.19) is regarded as a 



phenomenological prescription devoid of any first-principle motivation and cannot be implemented 



in the quantum theory as we did in Sec. IV C as it introduces a new length scale C. This modification 
is mandatory only for the situation in which, towards the singularity, fl x A x and fl p A p are of 
opposite signs, or equivalently dA/dr and dR/dr are of opposite signs according to (5. 15). 29 That 
is, at classical level, a given comoving cell is collapsing in the homogeneous {dQ,d^) directions 
while stretching in the radial (d x ) direction (or the other way around), reminiscent of the Kasner 
(vacuum) solution to the Bianchi I cosmology (see |17j for the analogy between the Bianchi I and 
the Kantowski-Sachs models). In the context of Bianchi I models, it has been shown in [35 1 that 
inclusion of (perfect fluid) matter with w < 1 has the effect of "isotropizing" the universe and loop 
quantum corrections take effect in the "isotropized phase" rather than in the "Kasner phase" if 
the matter content is abundant enough. 30 Likewise, it is anticipated that, in spherically symmetric 
models, with inclusion of abundant matter of w < 1, the solution is isotropized and loop quantum 
corrections take place in the isotropized phase. Therefore, for the cases with abundant matter 
content (collapsing black hole, Reissner-Nordstrom black hole, etc.), fl x A x and fi p A p approach the 
critical value ±7r/2 with the same sign and the original prescription in (4.25) can be used and the 
further modification in (5.18) other than I = may not be necessary. 



As commented in the end of Sec. VA the spherical symmetry is expected to break down in 
the vicinity of the quantum bounce from the perspective of the full theory of LQG. For the case 
devoid of or deficient in matter content, the quantum bounce takes place in the Kasner phase, in 
which the spacetime is highly anisotropic; consequently the breakdown could be exaggerated by 
the Kasner-like anisotropy and become even earlier and severer. Accordingly, the modification in 



(5.19), albeit an ad hoc prescription, can be interpreted as a correction which makes amends for the 



severe breakdown due to Kasner-like anisotropy, as C signals the length scale at which the weave 
state starts to become granular and / characterizes the degree of coherence of the weave state. 



29 Note that, towards a spacelike singularity, dA/dr and dR/dr dominate the terms involving spatial derivatives in 
(5.151, provided that the coordination and N x are well-behaved. Also see Footnote 28 

30 The analysis in [35] is based on the /i-scheme, instead of the /i'-scheme, but the conclusion that loop quantum 
corrections take effect in the isotropized phase when the matter content is abundant should still hold for the 
/i'-scheme. 
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C. Preliminary analysis of the quantum Schwarzschild spacetime 



At the heuristic level of effective dynamics, we treat the dyn amics as classical but governed 
by the holonomized Hamiltonian constraint given by (5.3) with the nth order holonomy 

corrections. The corresponding Hamilton's equations are given in (Dl). In the case of n = oo, we 
have 
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= fn{fJ> x A x ), 
= fn(fipA p ), 



(5.27a) 
(5.27b) 
(5.27c) 



by (5.2), (5.7) and (5.9). Taking (5.27) into (Dl), we obtain the Hamilton's equations when C is 
replaced by C^°°K 



In particular, (|2.45c|) and (|2.45d|), which in the classical case read as 
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(5.28a) 
(5.28b) 



are now replaced by 
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EP = Nf n (fL p A p ) ( S gn(E x ) 
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(5.29a) 
(5.29b) 



If we choose N x = 0, it is obvious that E x gets bounced as /n(Ax^x) flips signs, and E p gets 
bounced as f n (fi p A p ) flips signs. 

With the Hamilton's equations at hand, the numerical solution can be obtained by the method 
of finite- difference time-domain (FDTD) |36j. In the n = oo limit, however, the numerical method 
is hindered by the discontinuity appearing on the right hand sides of the Hamilton's equations. One 
can bypass this problem and still get faithful numerical solutions by keeping n large but finite. 31 



It is shown in [53] that, in the quantum theory of LQC, the evolution of coherent states follows smooth trajectories 
despite the kink of the connection variable with the n = oo holonomy corrections. In a sense, the discontinuity 
resulting from the n = oo limit is smeared by the quantum spreading. To reflect the smearing at the level of 
heuristic dynamics, one can simply truncate n to a large but finite value (the larger n effectively gives rise to the 
trajectory of a sharper coherent state). 
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Particularly, to obtain the loop quantum geometry of the maximally extended Schwarzschild 
spacetime, we choose the gauges 



N(t,x) 
N x (t,x) 



V32G 3 M 3 



IE" 



,1/4 



(5.30a) 
(5.30b) 



in agreement with (3.7), and accordingly use the classical solution (3.9) to yield the initial condition 



in the classical regime (say, for the t = spatial slice). That is, we assume that the weave state, after 
coarse-graining, manifests semiclassical traits (spreading, squeezing, etc.) in accordance with the 
Kruskal coordinates at large scale; this is a sensible prescription since the Kruskal coordinates give 
a well-behaved foliation which covers both the interior and exterior without introducing boundary 



terms, i.e. (2.53) satisfied with N± = (also see Footnote 28). With the gauges N and N x fixed and 



the initial condition given, the Hamilton's equations (Dl) in principle can be numerically solved 
for a given large n. The numerical task is however tremendously demanding on the numerical 
accuracy and thus it remains challenging and may require more sophisticated algorithms. 32 

Even though the numerical solution is extremely difficult, it still offers considerable insight 
about the natures of the loop quantum geometry of the Schwarzschild spacetime by inspecting the 



Hamilton's equation of the heuristic effective dynamics. Firstly, (5.29) implies that E x and E p get 
bounced as fn(fixA x ) and fn(fi p A p ) flip signs respectively. By the prescription (5.19), (5.26) further 



tells that E x and E p get bounced almost around the same epoch, provided that the parameter I is 
tuned to be close to but smaller than 3/4. This suggests that, for a well-behaved weave state which 
gives the extended Schwarzschild spacetime solution in the classical regime, the (black/white hole) 
singularity is resolved by the loop quantum corrections and replaced by a quantum bounce which 
bridges the classical solution to another classical phase. 

Secondly, the quantum effects are expected to become significant only for the region \E X \ < C 2 . 



However, (5.23) indicates that any spacetime point in the region of \E X \ > (2GM) 2 3> C 2 eventually 
receives quantum corrections in the late times (i.e. when t is large enough). In other words, the 
slices of fl x A x w ±7r/2 and ji p A p ~ = F 7r /2 inevitably intersect the event horizons. This seems 
to signal the breakdown (of the semiclassical treatment) of the quantum theory, as the spacetime 
curvature around the event horizon is fairly flat and should not incur any quantum corrections. This 
problem can be regarded as an indication that the inclusion of the Hawking radiation mechanism is 
mandatory. That is, only when the Hawking radiation is taken into account can the loop quantum 
theory be consistent. Since the Hawking radiation is not considered, the quantum corrections on 
the late-time horizons and beyond are not trustable. 

To summarize, assuming that the weave state is carefully chosen such that its semiclassical 
traits accord with the Kruskal coordinates and semiclassicality is upheld even across the putative 
singularity, two important consequences are observed: first, the classical singularity is resolved and 
replaced by a quantum bounce which bridges the black/white hole interior to a different classical 
phase; second, the Hawking radiation should be taken into account (on the late-time horizons) 
and thus the Hawking evaporation is expected. As both resolution of the classical singularity and 
inclusion of the Hawking radiation are expected, the complete quantum Schwarzschild spacetime 
is conjectured to have a global structure akin to that of the 2-dimensional black holes investigated 



32 The FDT D met hod is already very demanding on numerical accuracy even for solving the classical Hamilton's 
equations (2.451. For M « 10 2 — 10 3 mpi (mpi is the Planck mass), the numerical method requires 256-bit accuracy, 
which is about 77 digits, in order to approach the singularity close enough (until the quantum corrections take 
effect). The GNU Multiple Precision Arithmetic Library (GMP) [37] and the GNU MPFR Library [38] are needed 
for arbitrary-precision arithmetic and floating-point computations. To make sense of semiclassicality, M has to be 
much larger. The large M and higher order quantum corrections together demand even much higher accuracy. 
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in [28H30| . That is, the black hole is evaporated via the Hawking radiation, and the quantum 
spacetime is largely extended from the classical one as the classical singularity is resolved and the 
black hole interior is extended. The conjectured Penrose diagram is depicted in Fig.[3j Note that the 
black hole is simply evaporated; it is neither connected to a white hole nor to a baby black hole, 
as opposed to the results obtained in [13] (minisuperspace treatment) and [22] (midisuperspace 
treatment). As in the 2-dimensional quantum black holes in [28-30], the information that is 
classically lost in the process of Hawking evaporation is recovered, primary because the quantum 
spacetime is sufficiently larger than the classical counterpart. 

It should be noted that the scenario described above is expected only if the weave state manifests 
semiclassical traits in accordance with the Kruskal coordinates. Other scenarios are possible if one 
considers a weave state whose semiclassical traits are manifested differently (recall Footnote 28). A 
particular example is the solution obtained in |22j , which gives a singularity- free quantum spacetime 
akin to that in the fiQ- and /i-schemes of the minisuperspace (interior) treatment |11H13| and does 
not need to evoke the Hawking radiation. In [22], the 1-dimensional diffeomorphism gauge is 
partially fixed in such a way that the boundary conditions are given very similar to those in the 
Kruskal coordinates but one does not work exactly in Kruskal coordinates asymptotically. From the 
viewpoint of this paper, the gauge choice in |22] amounts to different semiclassical traits not exactly 
in accord with Kruskal coordinates and the free parameters in the solution could be interpreted as 
reflection of semiclassical traits (analogous to C and I in our treatment). 33 Finally, as commented 
in the end of Sec. |V A[ if semiclassicality breaks down severely in the quantum regime, the smooth 
weave state could simply descend into a granular spin network state. 



VI. SUMMARY AND DISCUSSION 



In Ashtekar's formalism for spherically symmetric midisuperspaces, the SU(2) internal gauge in 
the full theory is reduced to U(l) and the 3-dimensional diffeomorphism is reduced to 1-dimensional. 
By change of variables, the canonical structure in terms of (A x , A p ,n; E x , E p , P 11 ) decouples the 
U(l) degrees of freedom from the geometric/metric ones, and accordingly the Gauss constraint 
is completely decoupled from the diffeomorphism and Hamiltonian constraints. While (77; P 71 ) 
corresponds solely to the internal degrees of U(l) gauge, (E X ,E P ) gives the intrinsic geometry 
via (2.47) and (A X ,A P ), together with the lapse N and shift N x , gives the extrinsic geometry via 
(2.51). Particularly, the classical solution of the maximally extended Schwarzschild spacetime can 
be explicitly recast in terms of Ashtekar variables in accordance with the Kruskal coordinates. 

For the loop quantum theory of spherically symmetric midisuperspaces, we adopt a new ap- 
proach, which is different from that in the previous works |18H2(J] but directly in the same spirit 
of LQC, by considering cylindrical functions of the symmetry- reduced connection variables from 
the outset as in (4.1). Imposition of the Gauss constraint yields the [/(l)-invariant states as spin 
networks defined in (4.11). Further imposition of the diffeomorphism constraint via group averag- 
ing gives rise to spin- knots, which can be abstractly interpreted as ensembles of chunks of volumes 
arranged in a radial order with adjacent surfaces of areas without any reference to localization of 
the chunks and adjacent surfaces. Unlike the full theory of LQG, the volume and area eigenvalues 
are only partially quantized in the sense of Bohr compactification due to the fact that SU (2) is 
symmetry-reduced to U(l). 



Similarly, in the fio- and /Z-schemes of the minisuperspace (interior) treatment jllH13| . there appears to be an 
additional parameter in the effective solution which gives rise to the difference between the mass of the white hole 
and that of the black hole (essentially, the solution depends on the gauge choice the finite sized cell). This again 
can be interpreted as reflection of semiclassical traits of the given weave state (see |13| for more comments). 
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FIG. 3: (a): The Penrose diagram of the maximally extended Schwarzschild spacetime. (b): The (con- 
jectured) Penrose diagram of the loop quantum extension of (a). The putative singularities are resolved 
and replaced by the quantum bounces (shaded areas) and the generalized dynamical horizons (dashed lines) 
deviate from the classical event horizons of the black/white hole in the late/early time. The black and white 
holes are evaporated via the Hawking radiation (denoted by the arrows) . (More precisely, evaporation of the 
white hole is the time reversal of the black hole evaporation; the arrows pointing backwards in time denote 
forward radiation of antiparticles.) (c) and (d): Contours of constant R = \/\E x \ of (a) and (b) (smaller 
R in lighter color/shade). 

For the dynamics, we propose a new quantization scheme to construct a graph-preserving Hamil- 
tonian constraint operator, which changes only the coloring of spin network states and upon which 
one can impose the fundamental discreteness of LQG by hand a la the strategy of improved dynam- 
ics in LQC. The strategy of improved dynamics is demanded as for respecting the diffeomorphism 
invariance, which is absent but reduced to scaling invariance in the context of LQC. On the other 
hand, imposition of fundamental discreteness spoils the intimate matching between diffeomorphism 
and Hamiltonian constraints in the full theory of LQG, and as a consequence the Hamiltonian con- 
straint operators with different lapse functions do not commute, i.e. [C[iV], C[M]] ^ 0, and the 
resulting dynamics depends on foliation of the spacetime. However, this problem seems to be mit- 
igated at large scale, if one includes higher order holonomy corrections, which can be understood 
as a result of generic j representations for holonomies in the Hamiltonian constraint operator. 

While rigorous and complete construction of the quantum theory remains challenging, significant 
insights can still be obtained at the heuristic level of effective dynamics. Consistency of constraint 
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algebra regarding the Hamiltonian and diffeomorphism constraints fixes the improved quantization 
scheme to be of the form of (4.25), reminiscent of the /i'-scheme in LQC for the Kantowski-Sachs 
models. Moreover, consistency of constraint algebra regarding two Hamiltonian constraints (with 
different lapse functions) further demands the inclusion of higher order holonomy corrections as 
in (5.2) to the order of infinity and meanwhile fixes a ratio factor 2 in (4.25) for the improved 
quantization scheme. 

At the heuristic level of effective dynamics, it is suggested that the black hole singularity is 
resolved and replaced by a quantum bounce which bridges a classical solution to another classical 
phase, as the constraint algebra is exactly satisfied before and after the bounce if (4.25) and (5.2) 
with n = oo are adopted. As indicated in (5.13), however, the constraint algebra breaks down 
briefly during the transition period of the bounce, as fn{fixA x ) and fn(p> P A P ) flip signs at close but 
slightly different epochs. As a result, the two constraints C^ 00 -* = and D = are violated after 
the bounce, implying that the spacetime is no longer invariant under spatial diffeomorphism and 
change of spacetime foliation after the bounce. Nevertheless, if the characteristic parameter used 
to describe the solution in the classical regime is much larger than the Planck unit, the transition 
period (which is of the scale of Planck time) is so brief (compared to the characteristic time) 
that the violation is minuscule and therefore the heuristic effective dynamics still gives a reliable 
description of large-scale physics after the bounce. 

The fact that the constraints violate briefly during the transition period indicates that one 
cannot make sense of symmetry reduction by separating the degrees of freedom of the full theory 
into spherical (heavy) and non-spherical (light) ones in the vicinity of the quantum bounce. It 
suggests that, in the full theory of LQG, a given weave state which represents a smooth, spher- 
ically symmetric space at large scale will eventually manifest granularity of spin networks when 
approaching the quantum bounce. After the bounce, the weave state evolves to become smooth 
and spherically symmetric again provided that the transition period is very brief. In the case that 
fn(fixA x ) and fn(p, p A p ) do not flip signs at close epochs or even one of them never flips signs, 
however, a different scenario in the full theory is also possible. In this case, semiclassicality breaks 
down severely and the smooth weave state could simply descend into a granular spin network state, 
or other quantum effects (such as the Hawking radiation) should be taken into account to maintain 
the semiclassicality. 

At the heuristic level of semiclassical dynamics, further phenomenological modification can be 
prescribed as in (5.18) and (5.19). This modification is devoid of any first-principle motivation but 
phenomenologically can be interpreted as quantum corrections due to semiclassical traits (spread- 
ing, squeezing, etc.) of the weave state, as C signals the length scale at which the weave state starts 
to become granular and / characterizes the degree of coherence of the weave state. For the cases 
devoid of or deficient in matter content, the quantum bounce takes place in the Kasner phase and 
the phenomenological modification (5.19) is mandatory, if the semiclassicality is to be maintained 
across the putative singularity. On the other hand, for the cases abundant in matter content, the 
phenomenological modification (5.18) might be optional but the prescription with fine-tuned C 
and I could further mitigate the breakdown of semiclassicality. Semiclassicality of the weave state 
could be upheld even across the putative singularity, thus giving a bouncing scenario, if the weave 
state is carefully chosen (i.e. C and I are fine-tuned). 

The heuristic effective dynamics can be used to investigate the quantum extension of the 
Schwarzschild spacetime in accordance with the classical solution in the Kruskal coordinates. Al- 
though the full-fledged numerical method is extremely challenging, considerable insight can still 
be obtained by inspecting the effective Hamilton's equations: first, the classical singularity is re- 
solved and replaced by a quantum bounce which bridges the black/ white hole interior to a different 
classical phase; second, it is mandatory to take into account the Hawking radiation and thus the 
Hawking evaporation is expected. These lead us to conjecture that the complete quantum space- 
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time of the Schwarzschild solution resembles the quantum spacetime structure of the 2-dimensional 
black holes studied in |28H30j : the black hole is evaporated via the Hawking radiation, and the 
quantum spacetime is largely extended from the classical one as depicted in Fig. [3} The informa- 
tion that is classically lost through the Hawking evaporation is eventually recovered, because the 
quantum spacetime is sufficiently larger than the classical counterpart. 

It should be kept in mind, however, that other scenarios are still possible if the weave state 
exhibits different semiclassical traits. Particularly, the solution in |22j gives a complete quantum 
spacetime akin to that in the and /2-schemes of the midisuperspace treatment and does not evoke 
Hawking radiation. This very different scenario is not contradictory but in fact complementary 
to that conjectured in this paper, as the quantum evolution indeed depends on details (such as 
semiclassical traits) of the smooth weave state. The free parameters of the solution in [22J could be 
better understood, if one can rephrase the partial gauge fixing in [22] from the perspective of this 
paper. Furthermore, different scenarios (including those descending into granular spin networks) 
are in fact connected in the quantum theory in the sense that, at the quantum level, the outcome 
of a given weave state can branch into different scenarios with various probabilities. 

More insight about the natures of loop quantum corrections can be obtained if one applies the 
same formulation of this paper to other spherically symmetric models. For collapsing black holes, 
it would be very instructive to know how the formation of the trapped surface is altered as well as 
how the classical singularity is resolved by the quantum corrections. Moreover, for the Reissner- 
Nordstrdm (charged) black hole, one would have the chance to study the quantum corrections on 
the timelike singularity, which is not possible in the minisuperspace formalism. 
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Appendix A: Dimensions and density weights 

From the point of view of the 1-dimensional field theory over x, many variables transform as 
scalars or scalar densities under the 1-dimensional diffeomorphism on x. As noted in Footnote [6j 
it is very helpful to keep track of the density character and here we list the density weights for 
various variables in Table [TJ The transformation rule under diffeomorphism for scalar densities with 
different weights is given in Footnote [9j For convenience, we also list the dimensions for various 
constants and variables at the same time. We work in the units with c = 1 but keep G and h 
explicit. In this unit system, the Planck length is given by y/Gh/c 3 = VGh =: 4>i and the Planck 
mass by y/hc/G = y/h/G =: mp\. We keep the dimensions of the coordinates t and x arbitrary 
and denote the dimension of length as L. 



Appendix B: Notations for polar type variables 

In this paper, we adopt slightly different notations for polar type variables, which we believe 
are more succinct, self-explanatory and less misleading than those used in |18H22| . To reduce 
proliferation of notations, we also avoid defining new variables which are not used repeatedly. 



Furthermore, we also introduce the variables A r and P v as in Sec. II C the treatment of this 
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dimension 


density weight 


constants 




1 


— 




l\ x := Gh,A,£ 2 


L 2 


— 




mpi := \JhjG 


L/G 


— 


mass 


M 


L/G 


— 


Ashtekar variables 


A A 4 (n) 


1/x 


1 




A, A„ A A A^ 

■tt-ll A 1i sip, sip, Sip 


1 







a,/3,rj 


1 





triad variables 


E x 


L 2 







E 1 1 E 2 1 EP 


L 2 /x 


1 




pr) pi/ 


L 2 /x 


1 




q := ActEf = E X (EP) 2 


L 6 /x 2 


2 


spin network labels 




1 





discreteness variables 


e, fix 


X 


-1 




8, ftp 


1 









1 





3-metric components 


R 


L 







A 


L/x 


1 


constraints 


U 


h/x 


1 




D 


h/x 2 


2 




C, 


h/xL 


1 


£7(1) gauge angle 


X 


1 


(0) 


Lagrange multipliers 


UJ 3 ■ t 


lit 


(0) 




N x 


x/t 


(-1) 




N 


L/t 


(0) 


Hamiltonian density 


(w 3 • t) U+N X D+NC 


h/xt 


(1) 


Hamiltonian 


H 


h/t 





"The Lagrange multipliers as well as the U(l) gauge angle A transform as scalar densities with respect to the 
coordinate transformation x — > x(x) (together with t — > t = t) only in the sense that geometry of (SO(3, 1) bundle 
over) the spacetime is required to remain unchanged. As far as the diffeomorphism of the spatial slice alone is 
concerned, however, they are regarded as independent variables — i.e., they commute with the diffeomorphism 
constraint D. For A and variables involving Lagrange multipliers, we denote their density wights in parenthesis to 
indicate the subtlety. 

TABLE I: Dimensions and density weights of constants and variables. 



further change of variables is not explicitly carried out in |18H22j . The comparison of notations 
is given in Table |TT} Note that the spatial coordinates are denoted as (x,9,4>) in this paper, while 
they are denoted as (x, if) in [18^20] and (x, 9, (p) in [211 [22]; m the table, we denote them all as 
(x,9,4>). Also, the ditto mark " denotes that the case in the right column is identical to that in 
the middle column. 



Appendix C: Poisson bracket between Hamiltonian constraints 
1. Classical case 



At firs glance, it seems dauntingly tedious to derive (2.43f), as the Hamiltonian constraints 



iven in (2.42c) is rather complicated. A moment of reflection, however, tells that, essentially, only 
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in this paper 

(x,9,<f>) 


in [20] (and [HI [19]) 
{x,0 = O,<f> = <p) 


in [21] (and [22]) 

(x,9,<f> = <p) 




A X ,A 1 ,A 2 
E X ,E\E 2 


A X ,A X ,A 2 
E X ,E\E 2 


II 
II 




a,/3,ri 




II 




A P ,EP 


A^E* 


II 




Tl,T2,T 3 
Tl>,T 2 >,T 3 

Ty := T\ cos(a+/3) + t 2 sin(a+/3) 
ti> := -T\ sin(a+/3) + t 2 cos(a+/3) 

undefined 


Ai, A 2 , A 3 
A|,AtA 3 
A|; := Ai cos(a+/3) + A 2 sin(a+/3) 
K% := -Ai sin(a+/3) + A 2 cos(a+/3) 
A^ 1 := Ai cos /? + A 2 sin /3 
A^ := -Ai sin/3 + A 2 cos/3 


II 
II 

II 
II 


a 
b 


f}"^> pi p2 

-if =ri'/sind = -^ 


eaj, e^, sin 9 
— - -^Z 

L 4> ■— e x ~ 2E<t> 


p 3 9 J> 


c 


P p , P 13 defined in Footnote 11 
not used elsewhere 


P 4> .- ^E^ cos a 
PP := 2 A^E* sin a 


11 




A p := 2A p cos a 
P r > := 2A p EPsma 


2^K^ = 2A^ cos a 
pn = p0 ■- sin a 


:= 2^4^ cos a — 2^K^, 
n 




Ax := 4, + (a + /?)' = 7^ 


7^ = A, + (a + P)' 
undefined 


n 

undefined 





"In |18H21| . Ag and A^, are mistaken for each other. This typo is corrected here. 

''In [18H20| . A$ and A^ are mistaken for each other. This typo is corrected in Equation (36) of [21] and here. 
c Equation (29) of [21] reads as := — (e*)'/e^, which mistakes eg for as a typo. 



TABLE II: Comparison of notations for polar type variables. 



the Poisson brackets between Ashtekar connections and derivatives of conjugate triad variables will 
contribute. This is because anything like f dxdyf(x,y){A x (x),E x (y)} without derivatives inside 
the Poisson bracket is to be canceled out by the counter term J dxdyf(y,x){E x (x),A x (y)}. Thus, 
for example, we can compute the term 



A 2 p E p ^W\E x 'Ef" \ 

"4 7 2 vdW (EPf [y) j ~ A^^\EHx)\ (EP(y)) 



(Cl) 



272^^ (EP(y)f 
without worrying {A p (x), l/(E p (y)) 2 } 7^ 0. This largely simplifies the calculation and it turns out 



{C[N],C[M]} 



dxdyN(x)M(y){C(x),C(y)} 



dx 



sgxi{E x ) 
2G 7 {EP) 2 



j [A X E X E X ' - A p E x E p ' + A p E p E x '^j (MN' - NM') 
+A p E p E x (MN" - NM")} 



[dx (NM' - MN') ^^^ ( A'E P - A X E X ' 
2Gj J y J \E x \(E p ) 2 V p 



(C2a) 
(C2b) 



which leads to ( 2.43f ). In the abo ve, we have applied integration by parts on the term associated 
with A p E p E x ' in (|C2a[) to obtain \C2b\. 
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2. n — > oo case 



The same arguments used to derive (C2) also apply to the case 

of {cH[^],cW[M]}. That 

is, only the Poisson brackets between holonomized connections and derivatives of conjugate triad 



variables will contribute. Consequently, the only relations that matter are (p. 8a) and (5.8b). The 



factors fn(^xA x ) and fn(l^ p A p ) on the right hand sides of (5.8a) and (5.8b) will give rise to extra 



factors of the same kinds. For example, (CI) is modified accordingly as 



(A^f 



EP \/\E x \ E x 'E pl 



oo). 



EP(x) ^/\EHy)\E*'{y) 
^^/W(x)\ (EP(y)Y 

nr . f - 7 , A^\x)EP{x)^Wiy)\E x '{y) 

2G7/n M wjwm my)) 2 A ~ v) " 



^))\EP\y)) 



(C3) 



Taking care of all factors of fn(^xA x ) and f n (fJ, p A p ), we obtained the counterpart of (C2) for the 
case of n — > 00: 



{c^[N],C {oo) [M}} 



dx 



AG^iEPy i [V^ x A x )AWe*E*' - 2HU P A P )A^E*E<>' 

+ [fnffixAc) + Hfi p A p )) A p ^EPE x '] (MN' - NM') 
+2f n (fi x A x )A p ^EPE x (MN" - NM")}. 



Before and after the bounce, (C4) is largely simplified as 

1 



) ( A (00) ' pp _ A (oo) px' 

\E*\(EP) 2 \ p x 



, dx(NM' - MN') 
± : before/after bounce. 



(C4) 



(C5) 



By (5.12), we then have 



C^[N},C (oo) [M}} = ^ J dx (NM' - MN') ^ x ^ p)2 (D(x) - i/P") , (C6) 



which leads to (5.13). The symbol = is used in (C6) to indicate that the identity is valid only 



before and after the bounce but breaks down briefly during the transition of the bounce, because 
fn(fixA x ) and f n (fl p A p ) flip signs at close but slightly different epochs during the transition and 



thus the simplification in (C5) cannot apply. 



Appendix D: Hamilton's equations of heuristic effective dynamics 

At the level of heuristic effective dynami cs wi th the nth order holonomy co rrec tions, the corre- 
sponding Hamilton's equations is given by (|2.45|) with C replaced by in Kty. While (|2.45e|) 
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and (2.45f) remain unchanged, (2.45a)-(2.45d) are modified accordingly and take the forms as 

(Dla) 



A x = (2.45a) 



A x >A X \Ap >--Ap ^ 



-N sgn(E x ) 



7 



9^ 



+ sgn( J B x )- v 



dA 



(n) 



7 



A„ = (2.456) 



A x ^A x .An ^A 



(-0 



dE x ' 

(Dlb) 



-AT sgn(£ x ) 



iV 



SEP \2 7v ^ 



+ sgn(£ x )- v 1 



7 



= N\ sgn(E x 



EP = N\ sgn(E x 



\{n) 



dA, 



+ N x E a 



^W\At ] Af ] Ef \ 8A 



f(n) 



+ 



dEP ' 

(Die) 

(Did) 



(oo) 



For the case of n = oo, the Hamilton's equations above are large ly simplified: derivatives of A 
and Ap°°^ with respect to canonical variables are given by (5.27). Note that prior to the bounce, 
by (5.9) and (5.10), the Hamilton's equations in (Dl) for n = oo are identical to the classical 
counterparts in (2.45). 



With the prescription of (5.19) in particular, (5.27a) yields 



dAT ] 


k 


3E X 




dA^ 




dEP 








dE x 




p = 

dEP 



?r(4Z 2 -l) fc 



EP 



2^/A \E* 



n(2l-l) ( 
2^/A \E X 



before [i x A x reaches ±7r/2, 
r, after fi x A x , reaches ±7r/2, 

before p, x A x reaches ±7r/2, 



j72, after fi x A x , reaches ±7r/2, 



before /ipAp reaches =p-7r/2, 



(Ex) 



572, after /J, p A p , reaches =F7r/2, 



(D2a) 

(D2b) 

(D2c) 
(D2d) 



where we assume E x > for simplicity as E x does not flip signs throughout evolution. 



44 



[1] A. Ashtekar and J. Lewandowski, "Background independent quantum gravity: A Status report," Class. 

Quant. Grav. 21, R53 (2004) [gr-qc/0 404018 . 
[2] C. Rovelli, "Quantum Gravity," Cambridge University Press (2004). 

[3] T. Thiemann, "Modern Canonical Quantum General Relativity," Cambridge University Press (2007) 
[gr-qc/0110034| . 

[4] M. Bojowald, "Loop quantum cosmology," Living Rev. Relativity 11, 4 (2008); 

http://relativity.livingreviews.org/Articles/lrr-2008-4/. 
[5] A. Ashtekar and P. Singh, "Loop Quantum Cosmology: A Status Report," Class. Quant. Grav. 28, 

213001 (2011) |arXiv:1108.0893"l [gr-qc]]. 
[6] M. Bojowald, "Canonical Gravity and Applications: Cosmology, Black Holes, and Quantum Gravity," 

Cambridge University Press (2011). 
[7] D.- W. Chiou and M. Geiller, "Unimodular Loop Quantum Cosmology," Phys. Rev. D 82, 064012 (2010) 

[arXiv:1007. 0735 [gr-qc]]. 

[8] L. Modesto, "The Kantowski-Sachs space-time in loop quantum gravity," Int. J. Theor. Phys. 45, 2235 

(2006) [arXiv:gr-qc/0411032] . 
[9] A. Ashtekar and M. Bojowald, "Quantum geometry and the Schwarzschild singularity," Class. Quant. 
Grav. 23, 391 (2006) | arXiv:gr-qc/0509075] . 
[10] L. Modesto, "Loop quantum black hole," Class. Quant. Grav. 23, 5587 (2006) |! arXiv:gr-qc/0509078| . 
[11] L. Modesto, "Black hole interior from loop quantum gravity," Adv. High Energy Phys. 2008, 459290 
(2008) [gr-qc/0611043] . 

[12] C. G. Bdhmer and K. Vandersloot, "Loop quantum dynamics of the Schwarzschild interior," Phys. Rev. 

D 76, 104030 (2007) [arXiv:0709.2129~ | [gr-qc]]. 
[13] D.-W. Chiou, "Phenomenological loop quantum geometry of the Schwarzschild black hole," Phys. Rev. 

D 78, 064040 (2008) |arXiv:0807.0 665 [gr-qc]]. 
[14] A. Ashtekar, T. Pawlowski and P. Singh, "Quantum Nature of the Big Bang: Improved dynamics," 

Phys. Rev. D 74, 084003 (2006) [gr-qc/060703"9] . 
[15] D.-W. Chiou, "Effective dynamics, big bounces and scaling symmetry in Bianchi type I loop quantum 

cosmology," Phys. Rev. D 76, 124037 (2007) |arXiv:0710.0"416| [gr-qc]]. 
[16] A. Ashtekar and E. Wilson-Ewing, "Loop quantum cosmology of Bianchi I models," Phys. Rev. D 79, 

083535 (2009) |arXiv:0903.3397| [gr-qc]]. 
[17] D.-W. Chiou, "Phenomenological dynamics of loop quantum cosmology in Kantowski-Sachs spacetime," 

Phys. Rev. D 78, 044019 (2008) [arXiv:0803.3659] [gr-qc]]. 
[18] M. Bojowald, "Spherically symmetric quantum geometry: States and basic operators," Class. Quant. 

Grav. 21, 3733 (2004) [gr-qc/04070i7l . 
[19] M. Bojowald and R. Swiderski, "The volume operator in spherically symmetric quantum geometry," 

Class. Quant. Grav. 21, 4881 (2004) [gr-qc/04070i8j . 
[20] M. Bojowald and R. Swiderski, "Spherically symmetric quantum geometry: Hamiltonian constraint," 

Class. Quant. Grav. 23, 2129 (2006) |gr-qc/051H08] . 
[21] M. Campiglia, R. Gambini and J. Pullin, "Loop quantization of spherically symmetric midi- 

superspaces," Class. Quant. Grav. 24, 3649 (2007) j[gr-qc/0703135| . 
[22] R. Gambini and J. Pullin, "Black holes in loop quantum gravity: The Complete space-time," Phys. 

Rev. Lett. 101, 161301 (2008) |arXiv:0805.il87| [gr-qc]]. 
[23] K. V. Kuchaf, "Geometrodynamics of Schwarzschild black holes," Phys. Rev. D 50, 3961 (1994) |gr- 

qc/9403003. 

[24] D.-W. Chiou and L.-F. Li, "How loopy is the quantum bounce? A Heuristic analysis of higher order 
holonomy corrections in LQC," Phys. Rev. D 79, 063510 (2009) (arXi^090L1753 [gr-qc]]. 

[25] D.-W. Chiou and L.-F. Li, "Loop quantum cosmology with higher order holonomy corrections," Phys. 
Rev. D 80, 043512 (2009) |arXiv:0907.6640 [gr-qc]]. 

[26] M. D. Kruskal, "Maximal extension of Schwarzschild metric," Phys. Rev. 119, 1743 (1960). 

[27] G. Szekeres, "On the singularities of a Riemannian manifold," Pub. Mat. Debrecen 7, 285-301 (1960). 

[28] A. Ashtekar, V. Taveras and M. Varadarajan, "Information is not lost in the evaporation of 2- 
dimensional black holes," Phys. Rev. Lett. 100, 211302 (2008) |arXiv:0801.1811| [gr-qc]]. 

[29] A. Ashtekar, F. Pretorius and F. M. Ramazanoglu, "Surprises in the evaporation of 2-eimensional black 



45 



holes," Phys. Rev. Lett. 106, 161303 (2011) |arXiv:1011.6442| [gr-qc 



A. Ashtekar, F. Pretorius and F. M. Ramazanoglu, "Evaporation of 2-dimensional black holes," Phys. 
Rev. D 83, 044040 (2011) i arXiv:1012.0077| [gr-qc]]. 

A. Ashtekar, A. Corichi and P. Singh, "Robustness of key features of loop quantum cosmology," Phys. 
Rev. D 77, 024046 (2008) [ axXiv:0710.3565| [gr-qc]]. 

G. Date, "Absence of the Kasner singularity in the effective dynamics from loop quantum cosmology," 
Phys. Rev. D 71, 127502 (2005) [gr-qc/050500"2] . 

P. Singh and K. Vandersloot, "Semi-classical states, effective dynamics and classical emergence in loop 
quantum cosmology," Phys. Rev. D 72, 084004 (2005) [gr-qc/0507029] . 

M. Bojowald, "Large scale effective theory for cosmological bounces," Phys. Rev. D 75, 081301 (R) 
(2007) [gr-qc/0 608100 . 

D.-W. Chiou and K. Vandersloot, "The Behavior of nonlinear anisotropies in bouncing Bianchi I models 
of loop quantum cosmology," Phys. Rev. D 76, 084015 (2007) [arXiv:0707.2548| [gr-qc]]. 

B. Fornberg, "Generation of finite difference formulas on arbitrarily spaced grids," Mathematics of 
Computation, 51 (184): 699-706, (1988). 

The GNU Multiple Precision Arithmetic Library, http://gmplib.org. 
The GNU MPFR Library, |http://www.mpfr.org| 



46 



